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Introduction 


The  improvement  of  the  earth's  gravity  field  is  a 
continuous  task  of  geodesy.  In  the  past  years  one  procedure 
for  this  improvement  has  been  through  the  combination 
of  potential  coefficient  information  derived  through  the 
analysis  of  satellite  orbits  with  terrestrial  gravity  data 
(Kaula,  1966,  Rapp,  1978,  Lerch  et  als,  1979,  Gaposchkin, 
1980).  Many  of  the  solutions  that  were  made  were  carried 
out  with  the  gravity  data  defined  in  5°  equal  area  blocks. 
Such  data  was  to  be  the  dominant  source  of  information 
for  many  potential  coefficients  above  the  degree  of  satellite 
sensitivity.  Some  concern  was  expressed  in  Rapp  (1977a) 
about  the  information  (or  spectral)  content  within  such 
blocks.  In  that  study  it  was  pointed  out  that  information 
beyond  the  degree  irom  the  180°/8  rule  (0  is  the  block 
size  in  degrees)  may  be  in  the  data.  Since  there  was  some 
evidence  that  the  use  of  5°  anomaly  blocks  might  not  be 
the  best  size  to  use  when  higher  degree  solutions  were 
being  carried  out,  Rapp(1978)  carried  out  a  combination 
solution  with  l°x  1°  data.  This  solution  used  the  poten¬ 
tial  coefficients  and  their  accuracy  from  GEM9  (Lerch  et 
al  1979)  together  with  a  l°x  1°  anomaly  field  created  from 
the  merger  of  terrestrial  anomaly  data,  and  anomaly  data 
derived  from  the  Geos-3  satellite  altimeter.  The  combination 
solution  was  carried  out  with  a  priori  potential  coefficients 
taken  to  degree  12  only.  This  was  necessary  because  of 
the  large  amount  of  computer  time  that  would  have  been 
necessary  for  the  rigorous  adjustment  of  the  satellite 
data  and  the  6480C  l°x  1°  anomaly  set.  One  result  of  this 
adjustment  was  an  adjusted  set  of  64800  l°x  1°  anomalies 
that  were  then  converted  to  a  set  of  potential  coefficients 
complete  to  degree  180  using  the  usual  summation  formulas 
and  an  algoritum  described  by  Rizos  (1979).  This  set  of 
coefficients  was  described  by  Rapp  (1979)  and  used  for 
anomaly,  geoid  undulation,  and  deflection  of  the  vertical 
computation  as  described  by  Tscherning  and  Forsberg  (1981). 
The  results  in  this  latter  paper  (and  comparisons  with  a  180x 
180  potential  coefficient  by  Lerch  et  als  (1981))  indicated 
the  value  of  a  high  degree  and  order  potential  coefficient 
solution . 

One  of  the  unsatisfactory  aspects  of  the  Rapp  1978 
180x  180  solution  was  the  fact  that  satellite  implied  poten¬ 
tial  coefficients  above  degree  12  were  neglected  in  the 
solution.  This  implied  that  the  field  would  be  inadequate 
for  satellite  orbit  calculations.  We  then  examined  proce¬ 
dures  that  might  be  of  help  in  incorporating  the  higher 
degree  coefficients  without  significantly  changing  our 
computer  requirements.  The  results  are  described  in  Rapp 
(1980)  where  we  examined  combination  solutions  made  with 
the  rigorous  procedures  used  in  Rapp  (1978)  with  a  similar 


procedure  using  certain  assumptions  that  led  to  a  consid¬ 
erably  simplified  solution.  Computations  were  carried 
out  with  both  the  rigorous  and  approximate  approach  to 
generate  two  sets  of  potential  coefficients  to  degree  180. 

180.  One  solution  was  that  to  180  using  the  adjusted  anom¬ 
alies  described  in  Rapp  (1978)  based  on  a  combination  with 
the  GEli  9  coefficients.  The  two  field  were  then  compared 
to  find  the  average  percentage  difference  was  8.6°/,,  the 
root  mean  square  undulation  difference  was  ±80  c,  and 
the  root  mean  square  anomaly  difference  was  2  mgals.  The 
percentage  difference  was  6.0%  for  degrees  2  through  12, 
increasing  to  15.4%  for  25  through  36  and  then  descreasing 
to  about  8%  for  the  higher  degrees.  Considering  that  the 
data  noise  percentage  error  would  be  high  (>50%)  we  felt 
the  difference  between  the  rigorous  solution  and  the  approxi¬ 
mate  solution  was  small.  This  led  us  to  believe  it  would 
be  worth  while  to  pursue  an  approximate  combination  solution 
using  all  the  satellite  derived  potential  coefficient  infor¬ 
mation  as  well  as  the  latest  l°x  1°  data  that  might  be  avail¬ 
able. 

The  results  in  Rapp  (1980)  also  showed  the  problem 
of  using  5°  block  averages  in  the  combination  solutions. 

To  do  this  combination  solutions,  both  5°  and  1°  block 
data  wore  used  with  the  GEM  9  potential  coefficients.  We 
found  that  the  differences  between  the  coefficients  of 
the  two  solutions  increased  significantly  with  degree. 

For  example  at  degree  36  the  percentage  difference  reached 
74%  although  it  was  21%  at  degree  21.  These  differences 
are  caused  by  the  differences  in  the  spectral  content  of 
mean  values  of  different  sizes.  This  problem  was  examined 
by  Colombo  (1981)  who  carried  out  a  theoretical  and  numerical 
analysis  to  predict  the  sampling  error  caused  by  the  finite 
size  of  block  being  used.  The  results  obtained  by  Colombo 
agreed  quite  well  with  the  results  obtained  by  Rapp  (1980) 
in  a  purely  numerical  fashion. 

As  part  of  Colombo's  studies  it  became  apparent  that 
there  were  several  new  techniques  that  might  be  used  to 
obtain  an  optimal  combination  of  satellite  and  gravity 
data.  However  some  of  these  would  require  a  significant 
programming  effort.  However,  Colombo  had  indicated  procedures 
that  could  be  used  in  a  very  near  optimal  way  to  estimate 
potential  coefficients  to  a  high  degree  once  a  set  of  adjusted 
anomalies  on  a  global  basis  were  given. 

We  thus  decided  to  carry  out  a  new  combination  solu¬ 
tion  with  new  data,  and  using  the  approximate  type  of  com¬ 
bination  theory  that  would  incorporate  all  available  poten¬ 
tial  coefficient  information.  In  the  following  sections 
we  describe  the  theory  used,  the  data,  and  the  results. 
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The  Combination  Theory 


Let  C£m  ,  Sim  be  a  set  of  fully  normalized  poten¬ 
tial  coefficients  which  occur  in  the  following  description 
of  the  earth's  gravitational  potential  V  : 

CO  £ 

V(r,<?,X)  =  [1+  l  (~)2  l  (Czm  cosmX  +  S^m  sinmX)  • 

1=2  m=0 

•  Ptm(sin$) ]  (1) 


where: 
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are  the  geocentric  coordinates  of  a  point, 

is  a  scale  factor,  usually  taken  as  an  equatorial 
radius  , 

are  the  fully  normalized  potential  coefficients. 


If  we  are  given  a  set  of  anomalies 

da  ,  we  can  relate  the  coefficients  and  anomalies  with  the 
following  spherical  approximation: 
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where  y  is  an  average  value  of  gravity.  If  we  let 
N  =  Tr/e  where  0  is  the  mean  anomaly  block  size.  Colombo 
(1980,  p.71)  has  expressed  (2)  in  the  following  form: 
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Equation  (3)  is  referred  to  as  a  quadrature  formula  and 
U$,  is  called  a  de-smoothing  factor.  Colombo  showed  that 
the  quadrature  formula  can  give  results  for  the  estimation 
ot  the  potential  coefficients  almost  as  good  as  a  much 
more  complicated  optimal  estimation  procedure.  As  part 
of  the  use  of  the  quadrature  procedure  Colombo  recommended 
the  following  de-smoothing  factors: 


(2) 


(3) 


(4) 


if  0  t  £  £  N/3 
if  N/3  <  £  <  N 
if  £  >  N 


(5) 
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is  the  Pel linen/Meissl  smoothing  or  averaging  operator. 
We  have: 


=  cot  SU 


£ ( £+1) 


where  \i>a  is  the  radius  of  a  circular  cap  having  the  same 
area  as  Cjj  (Rapp,  1977a).  We  also  have: 

e£  =  1'— 5s~o  -2ltr  [P£-l(cOS^j)  -  P£+1(coS^)]  (7) 

where  P i  is  the  Legendre  polynomial  of  degree  £  .  A 
recursive  procedure  for  finding  £&  is  given  by  Sjoberg 
(1980). 

The  principle  of  the  combination  procedure  used  here 
was  first  discussed  by  Kaula  (1966).  The  method  is  based 
on  a  comparison  of  the  potential  coefficients  computed  from 
(2)  (or  (3))  with  those  values  derived  from  satellite  data 
with  an  adjustment  being  performed,  recognizing  all  the 
data  is  to  be  weighted,  to  obtain  a  consistent  set  of  potential 
coefficients  and  anomalies. 

We  briefly  describe  this  adjustment  process  as  follows: 

A  general  function  F  is  defined: 


F  =  F(L£a,  Lxa)  =  0  (8) 

where  are  the  adjusted  observations  and  Lx-  are 

the  adjusted  parameters.  A  linearized  observation  equation 
is  then  formed: 


Bt  Vz  +  Bx  Vx  +  W  =  0  (9) 

where 

B*  ‘  -iq  ■  B*  '  4i~  ■  w  ‘  F<LeV>  o°> 

where  are  the  actual  observations  and  Lx°  are  the 

observed  values  of  the  quantities  to  be  regarded  as  parameters 
(e.g.  the  potential  coefficients)  of  the  adjustment.  If 


and  Px  are  the  weight  matrices  for  the  observations 
and  parameters,  respectively,  we  have  for  the  correction 
to  the  observed  parameters,  Vx  : 


V  =  - C B '  M-1  Bv  +  Pv)-1  Bi  M_1  W 

X  X  XXX 


(11) 


with  the  corrections  to  the  observed  quantities  (e.g.  the 
gravity  anomalies),  " 

V£  -  P"1  B^  M’1  (Bx  Vx  +  W)  (12) 


where 


m  -  Ba  *r  Bi 


(13) 


In  our  case  we  have: 

F  -  K  -  Lx 


where  Ly  are  the  given  estimates  of  the  potential  coeffic¬ 
ients  (e.g.  the  GEM  9  coefficients)  and  L$  are  the  coeffic¬ 
ients  computed  from  (2)  or  (3)  with  the  observed  set  of 
gravity  anomalies.  In  this  case: 


Bx  =  I  (15) 

[B„  1  =  -fHI—  j  Y  ( <j> ,  X  )do  (16) 

1  JcJpc  y(£-1)  jo ,  .  nm 

*  J 

The  bracket  around  Br.  indicates  that  the  expression  on 
the  right  side  of  (16)  is  simply  one  element  in  the  Bn  matrix. 
We  note  that  (16)  applies  only  for  the  partial  derivations 
with  respect  to  potential  coefficients.  We  will  also  be 
interested  in  the  zero  and  first  degree  terms  of  the  spherical 
harmonic  coefficients  of  the  anomalies  themselves.  In  this 
case  (16)  would  be  written. 


'Vig--Ut  lol3 

Using  (15)  in  (8)  we  have: 
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Vx  =  -( ( B»  PT1  B')'1  +  px)_1  (B,  P71  B’ )_1  W 


(18) 


and  equation  (12)  reduces  to: 


VJ.  =  PV  B'l  P<  Vx  (19> 

From  equation  (18)  and  (19)  we  can  obtain  the  adjusted 
values  of  the  potential  coefficients  and  the  adjusted  anomalies: 


La  =  Lxo  +  Vx 

V  =  V  +  (20) 


Note  that  in  these  cases  we  have  taken  the  observed  values 
to  be  the  approximate  values  to  simplify  the  equation.  Having 
a  set  of  64800  l°x  1°  adjusted  anomalies  we  can  apply  (2) 
or  (3)  to  obtain  a  high  degree  spherical  harmonic  expression. 
The  resultant  coefficients  should  agree  exactly  with  the 
adjusted  potential  coefficients,  L*a  ,  with  the  higher 
degree  terms  representing  the  information  in  the  l°x  1°  anom¬ 
alies  . 

The  above  adjustment  equations  were  those  used  in  Rapp 
(1978).  The  computational  effort  in  evaluating  the  matrix 
B 2  P^1  B2  was  sufficiently  great  that  an  adjustment  to 
only  degree  12  was  made.  In  Rapp  (1980)  simplifications 
were  made  in  an  attempt  to  make  a  more  complete  solution. 

To  carry  out  this  simplification  we  assume  that  the  observa¬ 
tion  weights  for  the  anomalies  are  assigned  as  follows: 


rp  I  =  C.2.S(P 
L  Z1  mz 


(21) 


where  m  is  the  accuracy  of  a  l°x 1°  block  assumed  to  the 
same  for  all  blocks.  P2  is  thus  taken  as  a  diagonal  matrix. 
Assuming  that  Px  is  also  diagonal  (18)  becomes: 


[Vx] 


-fWl 

1  +  A[PX] 


where 


A 


n&  <f>  A  X  u  o 
(yU-1)}- 


(22) 


(23) 


Here  M  and 
of  the  block. 


AA  are  the  latitude  and  longitude  increments 
The  obvious  advantage  of  this  technique  is 
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that,  we  are  not  required  to  form  or  invert  the  normal  equation.- 
to  obtain  the  adjusted  coefficient  or  anomalies.  On  the 
other  hand,  to  gain  this  advantage'  we  must  assume  that  all 
the  l°x  1°  anomalies  have  the  same  accuracy  and  are  implicitly 
weighted  according  to  (21).  Examination  of  (22)  would  show 
that  we  are  really  just  computing  a  weighted  average  of 
the  a  priori  coefficient  and  that  implied  by  the  anomaly 
data . 


The  A  Priori  Potential  Coefficients  and  Their  Accuracy 


As  a  first  step  in  arranging  our  data  we  will  choose 
the  starting  potential  coefficients.  We  first  examine  the 
GEM  9  potential  coefficients  and  their  estimated  accuracy 
(Lerch  et  als  1979).  These  coefficients  are  based  on  satellite 
derived  information  only.  The  coefficients  are  complete 
to  degree  20  with  additional  coefficients  up  to  (30,  28). 

The  accuracy  of  the  geoid  undulation,  by  degree  for  the 
GEM  9  coefficients  is  given  in  Table  1.  The  overall  commis¬ 
sion  error  from  degrees  2  to  20  in  the  undulation  is  ±173  cm. 
The  geoid  undulation  error  by  degree  is  plotted  in  Figure  1  on 
page  8.  The  cumulative  error  is  shown  in  Figure  9  and  the 
percentage  error  in  Figure  11  on  page  33.  The  percentage  error 
is  defined  as  - - - 


*Et  - 


i  (ACjL+ASL) 

m=0 


(24) 


m=0 
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Since  the  GEM  9  potential  coefficient  set  is  now 
several  years  old  we  decided  to  estimate  a  second  set  of 
coefficients  on  the  basis  of  potential  coefficient  solutions 
that  now  seem  most  current.  In  doing  this  we  considered 
some  recent  solutions  that  were  derived  or  tailored  to  specific 
satellites  (such  as  Geos-3,  Seasat,  Lageos),  and  some  solutions 
that  derived  just  zonal  harmonic  coefficients,  or  resonance 
coefficients. 

The  only  separate  zonal  coefficients  that  were  considered 
were  the  odd  zonal  harmonics  that  were  estimated  by  King- 
Hele  et  al.  (1981b).  The  coefficients  (and  accuracies)  used 
from  this  paper  were  from  degree  3  to  degree  19. 

Coefficients  determined  by  resonance  analysis  were 
also  taken  from  several  sources.  Those  coefficients  for 
which  we  found  information  are  as  follows: 


1)  12th  order  terms  based  on  the  analysis  of  four  satellites 
(Reigber  and  Rummel,  1979).  The  recommended  solution  included 
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Table  1 


Accuracy,  By  Degree,  in  the  Geoid  Undulation 
Implied  by  Various  Geopotential  Models 


1  GEM  9  SET1  SET1  (adj) 


2 

3  cm 

2  cm 

2 

3 

12 

5 

5 

4 

9 

6 

6 

5 

20 

13 

13 

6 

16 

11 

11 

7 

29 

20 

18 

8 

25 

20 

17 

9 

38 

31 

22 

10 

35 

27 

19 

11 

49 

41 

23 

12 

43 

28 

18 

13 

50 

35 

20 

14 

45 

33 

18 

15 

51 

32 

17 

16 

47 

37 

18 

17 

53 

34 

17 

18 

55 

36 

16 

19 

53 

33 

16 

20 

53 

27 

15 

21 

26 

14 

22 

29 

15 

23 

26 

14 

24 

25 

14 

25 

28 

14 

26 

27 

14 

27 

22 

13 

28 

26 

13 

29 

22 

13 

30 

21 

12 

31 

21 

11 

32 

19 

11 

33 

19 

11 

34 

18 

10 

35 

21 

10 

36 

17 

10 
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all  12  order  coefficients  and  their  accuracy  from  degree 
12  to  degree  30. 

2)  13th  order  terms  from  Klosko  and  Wagner  (1981).  These 
terms  were  derived  from  the  analyses  of  13  satellites  yielding 
135  constraint  equations.  The  recommended  solution  yielded 
13th  order  coefficients  and  their  accuracy  from  degree  13 
through  30. 

3)  14th  order  terms  from  Reigber  and  Balmino  (1977),  and 
King-Hele  et  als  (1979).  The  Reigber  and  Balmino  solution 
used  data  from  9  satellites  to  derive  27  condition  equations 
which  led  to  a  recommended  set  of  coefficients  from  degree 

14  to  degree  30.  The  King-Hele  et  al  solution  estimated 
the  coefficients  from  degree  14  to  degree  22.  For  purposes 
of  merging  with  other  data  we  adopted  for  use  the  King- 
Hele  et  al  coefficients  to  degree  22,  and  the  Reigber  and 
Balmino  coefficients  from  there  to  degree  30.  The  coefficients 
were  also  checked  for  consistency  with  the  results  of  Kostelecky 
and  Klokocnik  (1979)  who  gave  odd  degree  coefficients  to 
degree  21 . 

4)  35th  order  terms  were  taken  from  King-Hele  and  Walker 
1981b).  In  these  computations  23  orbits  were  used  to  derive 
the  potential  coefficients  and  their  accuracy  from  degree 

15  through  degree  35. 

5)  30th  order  coefficients  were  taken  from  King-Hele  (private 
communication.  1981)  and  from  Kostelecky  and  Klokocnik  (1979). 
The  King-Hele  (30,  30)  coefficient  was  used  with  the  even 
degree  values  to  degree  42  from  Kostelecky  and  Klokocnik. 

We  should  note  the  care  in  which  the  use  of  the  coeffic¬ 
ients  derived  from  the  resonance  analyses  should  be  used. 

In  essence  a  lumped  coefficient  is  sensed,  and  the  coefficient 
is  separated  by  using  satellites  at  different  inclinations. 

If  sufficient  data  is  not  used  this  separation  will  not 

be  as  good  as  one  wants  or  needs.  In  addition  high  correlations 

may  exist  between  the  coefficients  of  a  given  solution. 

For  our  next  steps,  however,  we  will  assume  that  accuracy 
estimates  given  by  the  authors  are  valid  and  that  the  coef¬ 
ficient  estimates  are  independent. 

We  now  consider  more  general  solutions  that  may  provide 
information  towards  our  starting  potential  coefficient  set. 

The  ones  considered  were  as  follows: 

1)  the  GEM  9  set  which  was  previously  described; 

2)  two  solutions  that  were  developed  for  better  orbit  estima¬ 
tion  with  the  Seasat  satellite  (Lerch  et  als,  1981).  The 

two  solutions  examined  here  were  those  known  as  PGS-S2  and 
PGS-S4.  The  PGS-S2  solution  was  based  on  the  GEM  9  normal 
equations  plus  additional  laser  and  S-band  tracking  of  Seasat. 
The  S4  solution  added,  to  PGS-S2,  the  5°  anomaly  data,  Geos-3 
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and  Seasat  altimeter  data.  The  S-2  field  is  complete  to 
degree  30  with  additional  terms  to  degree  36,  while  S-4 
is  complete  to  degree  36.  We  will  choose  to  work  with  the 
S2  solution  because  of  concern  with  the  use  of  5°  block 
data  in  S-4  although  for  tracking  purposes  (with  Seasat) 

S-4  is  the  superior  solution; 

3)  the  PGS1331  field  (Marsh,  private  communication  1981) 
is  a  set  of  potential  coefficients  tailored  primarily  to 
the  orbit  of  the  Starlette  satellite.  This  was  started 
from  the  GEM  9  normal  equations  with  laser  observations 
on  Starlette  added  as  well  as  other  information  including 
satellite  altimeter  data,  and  some  laser  observations  on 
Lageous  .  The  PGS1331  field  is  complete  to  36,36  with  addi¬ 
tional  terms  to  48,43. 

4)  the  PGS  L-l  field  (Lerch  and  Klosko,  1981)  is  a  field 
built  through  the  addition  to  the  GEM  9  data  set  two  years 
of  laser  observations  on  Lageos.  This  field  is  complete 
to  degree  20  with  additional  terms  to  30,29.  The  L-l  field 
is  a  significant  improvement  over  the  GEM  9  field  at  the 
lower  (J,  £  5)  degrees.  A  set  of  standard  deviations  for 

each  coefficient  was  provided  by  Lerch  (private  communication, 
1981) . 

5)  the  Rapp  (1978)  180x180  field.  This  field  was  based 
on  a  combination  solution  with  GEM  9,  l°x  1°  terrestrial 
anomaly  data,  and  l°x  1°  anomalies  derived  from  the  Geos-3 
data  available  at  that  time. 

6)  the  NWL1G  solution  described  by  Anderle  (1979).  This 
field  was  one  tailored  for  Geos-3  orbits.  It  is  complete 

to  degree  13  with  additional  terms  to  28,27.  We  did  consider 
the  recent  potential  field  described  by  Gaposchkin  (1980) 
but  decided  not  to  use  it  because  essentially  the  same  data 
was  present  in  the  previously  described  solutions. 

Our  task  is  to  combine  the  potential  coefficient  infor¬ 
mation  available  in  the  above  field  to  obtain  a  starting 
set  of  potential  coefficients  and  their  accuracy  for  the 
combination  solutions.  Ideally  we  might  do  this  if  we  have 
the  variance-covariance  matrix  for  each  solution.  We  do 
not  have  such  information.  And  even  then  great  concern 
would  exist  when  combining  tailored  fields,  and  resonant 
coefficient  information. 

Before  merging  these  coefficient  sets  we  derived  a 
set  of  l°x  1°  anomalies  from  the  potential  coefficient  and 
compared  them  to  a  combined  terrestrial/Geos-3  l°x  1°  field. 

The  purpose  of  doing  this  was  to  see  if  one  or  more  field 
were  significantly  better  than  the  other  fields  in  comparisons 
with  anomaly  fields.  Such  comparisons  in  terms  of  the  mean 
square  anomaly  difference  are  given,  in  Table  2,  using  31210 
l°x  1°  "known"  values  having  an  estimated  standard  deviation 
smaller  than  11  mgals.  The  computations  were  done  for  two 
cases  of  the  maximum  degree  in  the  spherical  harmonic  expan¬ 
sion:  20,  and  36  (except  for  PGS-S2).  For  the  n  =  20  case, 

the  S-4  and  1331  fields  se®m  to  be  most  accurate.  For  the 
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case  of  n  =  36 ,  the  1331  field  seems  best.  We  should  note 
that  much  of  the  differences  seen  is  due  to  the  fact  that 
the  l°x  1°  data  has  much  high  frequency  information  not  nresent 
in  the  potential  coefficient  fields  being  tested.  We  conclude 
that  there  are  small  differences  in  the  solutions,  but  they 
are  perhaps  not  significant. 


Table  2 

Mean  Square  Difference  Between  l°x  1°  Anomalies 
Derived  from  a  Potential  Coefficient  Field, 
and  a  Terrestrial /Geos-3  Data  Set 


Field 


Mean  Square  Anomaly  Difference  (mgal2) 
NMAX=20  NMAX=36 


GEM  9 

365 

_ 

PGS  S-2 

362 

357* 

PGS  S-4 

353 

335 

PGS1331 

353 

331 

PGSL1 

365 

- 

"SET1' 

353 

334 

*to  degree  30  only 


We  therefore  choose  an  arbitrary  merge  procedure 
but  one  that  should  yield  realistic  estimates  of  the  coef¬ 
ficients  and  their  accuracy.  Specifically  we  found  a  coef¬ 
ficient  by  forming  a  weighted  average  of  the  PGS  S-2,  PGS1331, 
PGSL1 ,  and  the  miscellaneous  coefficient.  The  weighting 
for  the  LI  and  miscellaneous  coefficients  was  done  using 
the  standard  deviations  for  each  coefficient.  The  standard 
deviations  for  the  PGS1331  and  the  PGS  S-2  fields  were  taken 
to  be  0.9  that  of  the  corresponding  GEM  9  coefficient  where 
it  existed.  If  it  did  not  exist  we  used  the  standard  deviation 
of  the  miscellaneous  coefficient.  Since  the  1331,  S-2, 
and  L-l  fields  are  also  basically  dependent  on  the  GEM  9 
data  set,  so  then  does  our  final  result. 

To  compute  the  accuracy  estimate  of  a  particular  coef¬ 
ficient  we  choose  the  smallest  of  the  following  values: 

1)  the  standard  deviation  of  the  LI  solution; 

2)  the  root  mean  square  difference  between  the  weighted 
mean  value  and  all  other  corresponding  values; 

3)  the  standard  deviation  of  the  miscellaneous  coefficient. 

Again  the  above  procedure  is  somewhat  arbitrary.  It 
does  recognize  that  the  formation  of  a  weighted  mean  in 
the  procedure  used  here  may  not,  and  probably  does  not,  reduce 
the  given  error.  And  it  does  recognize  the  coefficients 
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that  have  good  agreement  between  the  various  data  sets  con¬ 
sidered.  Our  final  set  of  coefficient  formed  by  the  above 
procedures  will  be  called  "SET1". 

To  have  a  preliminary  check  on  the  starting  set  of 
coefficients  we  calculated  l°x  1°  anomalies  from  these  coef¬ 
ficients  and  compared  them  to  our  l°x  1°  terrestrial  data 
base.  The  mean  square  differences  are  shown  in  Table  2. 

We  see  that  there  is  no  significant  change  in  the  comparisons. 
The  geoid  undulation  error,  by  degree,  for  this  field  is 
given  in  Table  1  and  plotted  in  Figure  1.  The  overall  com¬ 
mission  error  to  degree  20  is  ±120  cm  as  compared  to  ±173  cm 
for  GEM  9.  The  commission  error  to  degree  36  is  ±152  cm.  The 
anomaly  degree  variances  have  also  been  computed  for  the 
GEM  9  and ( "SET1") coef f icients  and  the  results  given  in  Table 
3.  We  see  no  significant  difference  between  the  GEM  9  values 
and  those  implied  by  the  "SET1"  coefficients.  The  anomaly 
degree  variances  were  computed  as  follows: 

et  -  Y!u-1>'  J  (qm  +  5>m)  (25) 

m=u 


In  these  equation  the  even  degree  zonal  coefficients  are 
given  with  respect  to  the  reference  coefficients  implied  by 
an  ellipsoid  whose  flattening  is  1/298.257222.  A  more  accur¬ 
ate  formulation  of  degree  variance  computations  considering 
the  Bjerhammer  sphere  is  described  by  Jekeli  (1978).  A  plot  of 
the  anomaly  degree  variances  to  degree  180  for  the  adjusted 
field  is  given  in  Figure  12  on  page  34. 

Our  initial  procedure  was  to  carry  out  two  combination  sol- 
solutions:  one  with  the  GEM9  potential  coefficients  and 
the  other  with  the  SET1  coefficients.  Such  adjustments  were 
made  with  a  preliminary  l°x  1°  gravity  field.  Tests  with  the 
adjusted  fields  indicated  that  the  GEM9  combination  was  poorer 
than  the  combination  with  the  SET1  coefficients.  Therefore, 
we  decided  to  carry  out  only  one  combination  solution  and 
that  would  be  with  the  "SET1"  coefficients  and  their  accuracy. 
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Table  3 


Anomaly  Degree  Variances 

Implied 

By  Various 

Potential 

.  Coefficient 

Solutions  (mgal2) 

Si 

GEM  9 

"SET1" 

Adj .  with  "SET1" 

2 

7.56 

7.58 

7.58 

3 

33.66 

33.82 

33.84 

4 

19.63 

19.79 

19.81 

5 

20.87 

20.85 

20.70 

6 

19.04 

19.42 

19.33 

7 

19.45 

19.60 

20.10 

8 

11.73 

11.30 

11.02 

9 

11.50 

11.37 

11.09 

10 

10.1 

10.0 

9.77 

11 

6.8 

6.4 

6.85 

12 

3.7 

3.2 

3.2 

13 

6 . 6 

7.0 

7.2 

14 

4.0 

3.4 

3.0 

15 

3.3 

3.0 

3.0 

16 

2.3 

2.9 

4.8 

17 

2.1 

2.5 

3.5 

18 

3.3 

3.1 

3.3 

19 

3.0 

2.8 

3.0 

20 

2.3 

1.8 

1.9 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 


1.7 

2.7 

1.8 
1.7 
1.6 
1.4 

1.4 

2.1 

1.5 

1.6 
1.2 
0.9 
1.3 
1.0 
1.7 
0.7 


2.4 

3.7 

2.4 
2.1 
2.6 

1.7 
1.6 
2.3 
1.9 

2.5 

1.6 
1.6 
2.0 

2.8 
2.3 
1.7 


The  A  Priori  Gravity  Anomalies  and  Their  Accuracy 


The  l°x  1°  gravity  anomalies  to  be  used  in  the  combin¬ 
ation  solution  were  based  on  a  merging  of  a  terrestrial  data 
file  and  a  set  of  anomalies  derived  from  Seasat  anomaly  data. 

The  terrestrial  data  used  was  based  on  a  preliminary 
update  of  a  tape  known  as  the  October  1979  tape.  This  update 
used  ten  new  data  sources.  Of  special  importance  was  the 
revision  of  approximately  400  anomalies  in  the  southern  part 
of  Africa.  The  updated  tape  contained  42585  anomaly  values. 

The  location  of  these  anomalies  is  shown  in  Figure  2. 

The  Seasat  altimeter  data  was  adjusted  by  Rowlands 
(private  communication,  1981)  to  the  point  where  the  average 
crossover  discrepency  was  ±28  cm.  This  data  was  then  used 
to  estimate  37905  l°x  1°  anomalies  and  undulations  some  of 
which  were  on  land  and  thus  not  reliable.  The  location  of 
this  data  is  shown  in  Figure  3. 

The  two  anomaly  sets  were  then  merged  together  to  form 
a  set  of  56761  anomalies  which  are  shown  in  Figure  4.  In 
this  merger  the  following  criteria  were  used: 

1.  For  land  blocks,  use  only  the  terrestrial  data; 

2.  For  oceanic  blocks,  use  the  terrestrial  data  if  the  standard 
deviation  was  i±5  mgals; 

3.  For  an  oceanic  block,  bordering  land  use  terrestrial  data 
when  it  exists'; 

4.  Use  terrestrial  data  (when  it  exists)  in  the  Mediterranean 
Sea  area  due  to  possible  tide  problems  in  the  Seasat  analysis. 

5.  For  all  other  blocks  use  the  Seasat  derived  anomalies 
when  it  exists.  The  location  where  the  Seasat  anomalies 
(33905  values)  are  located  is  shown  in  Figure  5. 

All  anomalies  were  referred  to  the  Geodetic  Reference 
System  1980  (Moritz,  1980).  Before  any  altimeter  anomalies 
were  used  in  the  merger  the  atmospheric  correction  of  -0.87 
mglas  was  applied  to  the  anomaly  predicted  using  least  squares 
collocation . 

This  new  data  set  was  compared  to  our  earlier  merger 
with  Geos-3  data  (Rapp,  1980).  For  the  52972  common  values 
wer  found  a  root  mean  square  difference  of  ±7.5  mglas,  with 
a  maximum  difference  of  127  mgals.  The  root  mean  square 
standard  deviation  of  all  anomalies  on  the  new  merged  tape 
is  ±10.0  mgals. 

In  the  application  of  equation  (2)  or  (3)  to  the  anomaly 
data  a  global  estimate  of  64800  l'x  1°  anomalies  is  needed. 

The  anomaly  values  remaining  from  the  56751  "known"  values 
were  set  to  the  values  implied  by  the  "SET1"  potential  coef¬ 
ficients  to  degree  36.  If  one  considers  the  standard  deviation 
of  these  values  to  be  ±30  mgals,  the  root  mean  square  standard 
deviation  of  all  64800  l°x  1°  anomalies  is  about  ±15  mgals. 

For  purposes  of  evaluating  m  in  equation  21  we  choose  the 
average  standard  deviation  to  be  ±20  mgals  to  allow  for  possible 
anomaly  errors  in  the  several  largely  unsurveyed  areas. 
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The  Atmospheric  Correction  Terms  for  the  Anomalies 


The  gravity  anomalies  described  above  have  not  been 
corrected  for  the  effect  of  the  mass  of  the  atmosphere. 

In  order  to  properly  use  equation  (2)  or  (3)  such  corrections 
must  be  applied (Moritz,  (1974),  Rummel  and  Rapp  (1976), 

Rapp  (1977b)). 

If  Ag  are  the  gravity  anomalies  referred  to  an  ellip¬ 
soid  where  the  total  mass  of  the  atmosphere  has  been  condensed 
on  its  surface,  the  proper  anomaly  to  use  in  (2)  or  (3) 
is: 


Ag°  =  Ag  +  6gA 


(26) 


where  6gA  is  the  atmospheric  correction  which  has  been 
tabulated,  as  a  function  of  elevation,  most  recently  in 
Moritz  (1980a) .  (The  sign  convention  used  in  (26)  is  opposite 
to  that  used  in  earlier  papers  but  is  consistent  with  that 
used  in  Moritz  (1980b)).  SgA  is  0.87  mgals  for  ocean  blocks 
and  land  blocks  whose  elevation  is  zero.  For  blocks 
whose  mean  elevation  is  5000  m.  the  correction  is  0.47  mgals. 
Although  this  correction  is  small  in  terms  of  potential 
coefficient  effects,  it  is  easy  to  make  and  can  remove  a 
source  o^  systematic  error.  Such  corrections  were  made 
for  our  final  data  set. 


The  Ellipsoidal  Correction  Problem 


Equations  (2)  or  (3)  represent  spherical  approximations 
to  a  more  accurate  solution  considering  the  reference  surface 
to  be  ellipsoidal.  Tests  described  by  Rapp  (1977a)  using 
correction  terms  from  equations  of  Lelgemann  (1973)  showed 
that  the  corrections  were  small  and  could  be  neglected. 
Recently  Pellinen  (1981)  has  re-examined  the  Lelgemann  solu¬ 
tion  and  derived  new  corrected  formulas  that  are  in  the 
same  form  as  Lelgeman's.  We  can  write  from  Pellinen  (1981): 


where: 
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(  32, 2 -114  +14  )  (  Z  -m-1 )  ( 4  -pi ) 

2(£-l)(  24 -  3)(  22.-1,/ 

2f +  22jn:  -22. 3  -42nr  +9g. 2  +8m:  +9?  -8 

2(4-1)  (2?. +3  )(2M  )  (28) 

L2+52  +  2)(2+m+2)( 2+m+l) 

2(4-l)(24+5)(24+3) 


ucic  u^£ra  »  are  t^ie  correcti°ns  to  the  coefficients 

that  woula  be  computed  from  (2)  if  geodetic  latitudes  were 
used  in  the  computations.  These  correction  terms  were 
evaluated  for  the  GEM  10B  potential  coefficients  which 
are  complete  to  degree  36.  The  maximum  correction  we  found 
was  0.015 x 10-6  for  C2,2  •  Other  terms  were  considerably 
smaller.  For  example,  at  degree  22,  the  root  mean  square 
correction  was  ±0.0009 x  10-6  .  This  value  is  consistent 
with  the  expected  value  predicted  from  3xlO”8/4  given 
by  Pellinen  (1981)  which  was  derived  using  the  Kaula  10_5/22 
decay  rule  for  fully  normalized  potential  coefficients. 

We  have  decided  not  to  use  these  correction  terms. 
Instead  we  will  use  geocentric  latitudes  in  all  computations 
involving  (2)  or  (3).  This  does  not  eliminate  all  concerns 
about  the  ellipsoidal  and  spherical  reference  surfaces. 
However,  at  low  degrees  where  this  effect  will  be  greatest, 
the  adjusted  coefficients  will  be  dominately  determined 
by  the  satellite  determined,  a  priori  coefficients.  At 
the  higher  degrees  the  effect  of  the  spherical  approximation 
is  at  least  an  order  of  magnitude  below  the  uncertainty 
in  our  coefficients. 
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Results  of  the  Combination  Solution-General 


The  combination  solution  was  made  with  SET1  a  priori 
coefficients  and  the  64800  l°x  1°  starting  anomalies.  The 
first  result  of  the  adjustment  were  the  adjusted  l°x  1°  anomalies 
and  the  adjusted  potential  coefficients  corresponding  to 
the  a  priori  values.  A  number  of  comparisons  and  computations 
can  be  made  for  this  adjusted  data. 

We  consider  first  the  accuracy  of  the  geoid  undulation, 
by  degree,  implied  by  the  accuracy  of  the  adjusted  coefficients. 
The  values  were  given  in  Table  1  on  page  9  and  are  plotted 
in  Figure  1  on  page  8.  The  overall  error  from,  degree  2 
through  36  was  ±152  cm  for  the  starting  coefficients  and 
±87  cm  after  the  adjustment. 
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We  have  also  compared  the  adjusted  coefficients  with 
the  starting  coefficients  to  see  the  magnitude  and  location 
(by  degree)  of  the  changes.  These  comparisons  are  in  Table  4. 


Table  4 

Differences  Between  the  A  Priori  Coefficients 
and  the  Adjusted  Coefficients 


Solution  with  "SET1" 


l 

X 

AN (cm) 

6Ag(mgals) 

2 

0 

0 

.0 

3 

0 

1 

.0 

4 

0 

2 

.0 

5 

1 

5 

.0 

6 

1 

8 

.1 

7 

6 

27 

.2 

8 

9 

26 

.3 

9 

15 

40 

.5 

10 

14 

32 

.4 

11 

30 

52 

.8 

12 

19 

20 

0.3 

13 

24 

36 

0.7 

14 

52 

45 

0.9 

15 

35 

28 

0.6 

16 

42 

40 

0.9 

17 

42 

32 

0.8 

18 

39 

27 

0.7 

19 

47 

29 

0.8 

20 

55 

26 

0.8 

21 

47 

24 

0.7 

22 

41 

25 

0.8 

23 

47 

22 

0.7 

24 

47 

19 

0.7 

25 

53 

23 

0.9 

26 

57 

19 

0.8 

27 

59 

19 

0.8 

28 

49 

18 

0.7 

29 

76 

24 

1.0 

30 

63 

22 

1.0 

31 

81 

23 

1.0 

32 

65 

17 

0.8 

33 

62 

18 

0.9 

34 

61 

20 

1.0 

35 

61 

18 

0.9 

36 

79 

19 

1.0 

In  this  table  we  have  given  the  average  percentage  differ¬ 
ence,  the  root  mean  square  undulation  and  anomaly  difference 
by  degree.  We  see  the  percentage  changes  are  greatest  at 
the  high  degrees.  This  is  primarily  caused  by  a  higher  a 
priori  standard  deviation  at  such  degrees. 

The  anomaly  degree  variances  for  the  adjusted  solutions 
are  given  in  Table  3  on  page  14.  The  most  significant  changes 
occur  at  the  higher  degrees  with  the  new  solution  having 
somewhat  larger  values  than  the  starting  values. 

The  l°x  1°  adjusted  anomalies  were  compared  to  the  adjusted 
anomalies  of  the  earlier  (Rapp,  1978)  solution.  The  mean 
difference  (new-old)  was  0.5  mgals,  the  root  mean  square 
difference  was  ±11  mgals  and  the  maximum  difference  of  215  mgals 

We  also  examined  the  magnitude  of  the  residuals  with 
respect  to  the  standard  deviations  of  the  64800  l°x  1°  anomalies, 
although  such  standard  deviations  were  not  used  in  the  solutions 
The  results  are  given  in  Table  5. 


Table  5 


Root  Mean  Square  Residuals  as  a  Function 
of  the  Initial  Anomaly  Standard  Deviation 


Range  of  Initial 
Anomaly  Std.  Dev. 
(mgals) 

Adjustment  with 

SET1  (m  =  ±20  mgals) 

1-5 

±3.0  mgals 

6-10 

2.7 

11-15 

4.5 

16-20 

4.9 

21-25 

5.4  V 

26-30 

3.5 

31-35 

5.0 

We  see  from  Table  5  a  clear  tendency  for  the  residuals 
to  be  larger  where  the  anomaly  standard  deviations  are  larger. 

In  Table  6  we  are  given  the  number  of  residuals  having 
a  given  absolute  magnitude  value. 


Table  6 


Number  of  l°x  1°  Residuals 
Having  a  Specified  Range 


Range 

(mgals) 

Adjustment 
with  "SET1" 

0  to 

2 

30884 

2  to 

4 

19855 

4  to 

6 

8160 

6  to 

8 

3363 

8  to 

10 

1318 

10  to 

12 

664 

12  to 

14 

289 

14  to 

16 

116 

16  to 

18 

84 

The  location  of  the  3827  residuals  whose  magnitude  is 
greater  than  7  mgals  is  shown  in  Figure  6.  The  large  residuals 
appear  to  be  generally  (but  not  always)  correlated  with  the 
locations  where  the  l°x  1°  anomalies  have  been  geophysically 
predicted.  These  latter  values  (6413  on  our  latest  terrestrial 
tape)  are  shown  in  Figure  7. 


The  adjusted  anomalies  from  each  solution  were  developed 
into  spherical  harmonic  coefficients  using  subroutine  HARMIN 
(Colombo,  1981)  and  the  optimum  quadrature  weights  given 
in  equation  (5).  The  expansions  were  made  to  degree  300 
but  results  only  to  degree  180  will  be  discussed  here.  This 
is  because  the  quadrature  weights  from  (5)  have  a  sharp 
discontinuity  at  degree  180.  For  example,  at  degree  180 
( =180° /1° ) ,  the  value  of  3^  is  0.65  so  that  m s o  is 
0.65  but  n i a i  is  1.  An  improvement  is  needed  in  the  deter¬ 
mination  of  the  optimum  quadrature  weights  beyond  the  Nyquist 
frequency. 

The  spectrum  of  the  potential  from  the  expansion  of 
the  adjusted  anomalies  is  shown  in  Figure  8.  Here  the  unitless 
spectrum  is  compute  from: 


V*(AU) 


CO 


l 

m=0 


’!m' 


(29) 


An  analysis  of  this  spectrum  and  isostatic  compensation 
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Figure  8 


POTENTIAL  DEGREE  VARIANCES 

9 

v'(  U)  -  «=*  +  S2m) 


fROM  Kaula's  10"5/*2  Rule 


(Rapp,  1979) 


is  given  in  Rapp  (1982).  For  comparison  purposes  we  have 
also  plotted  the  spectrum  implied  by  the  Kaula  rule  of  thumb, 
and  the  spectrum  implied  by  an  earlier  180  model  (Rapp,  1979). 


Potential  Coefficient  Comparison 


The  potential  coefficients  of  this  combination  solution 
have  been  compared  to  a  number  of  other  coefficient  sets.. 
These  results  are  given  in  Table  7. 


Table  7 

Potential  Coefficient  Comparisons 
with  the  Adjusted  Coefficients  of  this  Report 


Average 

RMS  Coeff.  Percent  RMS  Und.  RMS  Anomaly 


Solution 

NMAX 

Difference 

Diff . 

Diff. 

Difference 

GEM9 

20 

±13.5 

X 

10"9 

34 

±180  cm 

±  3.6 

mgals 

GEM 1 0B 

36 

±  6.9 

X 

lO"9 

47 

±162  cm 

±  5.0 

mgals 

GRIM3 

36 

±13.6 

X 

10*9 

60 

±319  cm 

±  6.5 

mgals 

Rapp  1979 

180 

±  1.6 

X 

10~9 

42 

+186  cm 

±  9.1 

mgals 

GEM10C 

180 

±  1.7 

X 

10"9 

81 

±201  cm 

±17.3 

mgals 

The  percentage  difference  between  the  Rapp  1979  field 
and  the  newer  field  is  10%  at  degree  7,  23%  at  degree  10, 
increasing  slowly  to  a  60%  difference  near  degree  180. 

The  percentage  difference  between  the  GEM10C  field  (Lereh 
et  als.,  1981)  and  the  new  field  is  7%  at  degree  7,  14%  at 
degree  10,  increasing  slowly  reaching  120%  near  degree  180. 

Of  special  interest  is  the  very  large  anomaly  difference 
between  the  two  solutions  of  ±17  mgals.  Much  of  this  can 
be  related  to  the  combination  techniques  used  and  the  basic 
l°x  1°  data  sets  which  is  somewhat  different  for  each  solution. 


Adjusted  Coefficient  Accuracy 


The  standard  deviation  of  the  adjusted  coefficients 
is  part  of  the  adjustment  solution.  When  carrying  out  the 
high  degree  expansions  it  is  necessary  to  assign  approximate 
accuracy  estimates  to  the  coefficients  that  are  not  part  of 
the  adjustment.  To  do  this  we  consider  two  error  components: 
the  first  is  due  to  the  data  noise,  while  the  second  is  the 
sampling  error  (Colombo,  1981)  due  to  the  finite  size  of 
the  l8x  1°  blocks . 
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The  standard  deviation  of  a  fully  normalized  potential 
coefficient  of  degree  l  based  on  anomaly  data  in  a  block 
of  size  8 (radians)  can  be  written  as  (Rapp,  1972): 


m(C,S). 


'NOISE 


_  m(A  g)  0 

/if 


For  our  case  0°  =  1°  and  m(Ag)  is  20  mgals. 

The  sampling  error  has  been  modeled  by  Jekeli  based 
on  computations  of  Colombo.  This  error  can  be  expressed 
as  Colombo,  (1981,  eq.  3.10): 


m(C,S)jl  =  -gjg-  ([(-16.19570  (-J-) 

SAMPLE 

+  30.34506  (-J-)  +  40.29588]  (-|-)2) 


where  is  the  root  mean  square  coefficient  of  degree 

l  ,  and  N  is  the  Nyquist  degree  which  is  180  for  the 
l°x  1°  data,  Equation  (31)  was  evaluated  by  first  computing 
a i  based  on  the  given  coefficients.  The  two  standard  devia¬ 
tions  were  then  quadratically  added  to  yield  a  composite 
standard  deviations  for  those  coefficients  not  part  of  the 
actual  adjustment.  This  merged  set  of  coefficients  were 
then  used  to  compute  several  accuracy  estimates  of  interest. 
Examples  of  the  percentage  error  for  the  data  noise  and 
the  sampling  error  are  given  in  Table  8. 


Table  8 

Coefficient  Percentage  Error  Due  to  Anomaly  Data 
Noise  (  ±20mgals)  and  the  Sampling  Error  ( l°x  1°  Blocks) 


Degree 

Percent 
Data  Noise 

Error 

Sampling  Error  j 

50 

52% 

n  i 

75 

72 

9  I 

100 

86 

16  1 

125 

98 

26 

150 

127 

38  ' 

175 

138 

52  j 

We  see  that  the  sampling  error  is  considerably  smaller 
than  the  data  noise  at  the  lower  degrees  but  increases 
at  the  high  degrees.  In  Figure  9  we  have 
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plotted  the  geoid  undulation  error  cummulativley  for  the  GEM9 
coefficient  set  by  degree  and  cumulatively  for  the  adjusted 
coefficients  of  this  report. 


In  Figure  10  we  have  repeated  the  information  of  Figure  9 
but  for  anomalies  that  would  be  calculated  from  the  spherical 
harmonic  expansions.  Also  shown  is  the  cummulative  anomaly 
error  for  the  a  priori  coefficient.  We  next  show  in  Figure 
11  the  percentage  error  several  solutions.  We  see  that  the  per¬ 
centage  error  of  the  adjusted  coefficients  reaches  100%  near 
degree  120.  The  increasing  percentage  error  beyond  this 
reflects  the  fact  that  the  coefficients  decay  faster  than 
the  error  in  the  coefficient  determinations.  Although  the 
percentage  errors  are  large  at  high  degrees,  the  coefficient 
set  as  a  whole  contains  significant  information. 


Anomaly  Degree  Variance  Models 


In  a  previous  section  (see  Table  3,  p.  14)  we  have  dis¬ 
cussed  anomaly  degree  variances  from  various  solutions  to 
degree  36.  It  is  now  appropriate  to  consider  these  anomaly 
degree  variances  for  the  adjusted  coefficients,  of  GEM10C 
and  GRIM3  (Reigber,  Balmino,  Moynot ,  and  Muller,  1981)  which 
are  complete  to  degree  36.  These  degree  variances  are  plotted 
in  Figure  12  along  with  that  implied  by  Kaula's  rule.  The 
values  from  the  GRIM3  model  are  higher  after  degree  20  than 
the  values  from  GEM10C,  and  the  new  model.  This  is  also 
reflected  in  the  summation  of  the  anomaly  degree  variances 
from  degrees  2  to  36.  This  summation  is  as  follows:  GEM10C, 
239  mgals;  GRIM3,  265  mgals;  new,  228  mgals.  The  anomaly 
degree  variances  for  the  GEM10C  field  are  lower  at  the  higher 
degrees  than  those  implied  by  the  new  field.  This  may  be 
due  to  the  analysis  technique  and/or  the  data  used. 


Orbit  Calculations  with  the  New  Models 


The  new  potential  coefficient  set  has  been  tested  in 
a  number  of  different  satellite  dependent  calculations.  This 
testing  has  been  to  see  how  badly  the  new  model  does  in  orbit 
work.  This  pessimism  stems  from  the  use  of  a  coefficient 
set  that  has  not  been  optimized  for  a  certain  satellite. 

We  have  tried  to  generate  a  general  field  that  would  be  useful 
for  both  terrestrial  applications  and  satellite  applications. 

The  tests  to  be  described  were  made  at  NASA's  Goddard 
Space  Flight  Center  by  Frank  Lerch  and  James  Marsh.  The 
tests  are  far  from  exhaustive. 

The  first  test  was  made  with  Lageos  in  a  25  day  arc. 

The  computations  were  made  with  the  complete  25  day  arc, 
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FIGURE  12 


Anomaly  Degree  Variances  For  GRIM  3,  GEM10C, 
Kaula's  Rule  And  For  The  Solution  Of  This  Report 


— Kaula's  Rule 

.GRIM  3 


Jhis  Report 


and  with  5  days  deleted  out  of  the  middle  of  the  arc.  Then 
the  root  mean  square  orbit  position  difference  and  the  root 
mean  square  residual  for  data  within  the  5  day  segment  were 
computed  using  several  potential  coefficient  models.  These 
results  are  given  in  Table  9. 


Table  9 

Computations  Using  Lageos  with 
Different  Potential  Coefficient  Models 


GEM 9  PGSL1  Rapp  (B) 

RMS  Orbit  ±1.92  m  ±0.81  m  11.06  m 

Post.  Diff. 

RMS  Obs.  ±22  cm  ±19  cm  ±20  cm 

Residual 

with  5  day 

gap 


The  next  test  involved  the  computation  of  10  baselines 
involving  5  laser  stations  and  two  different  sets  of  laser 
observations.  The  root  mean  square  differences  between  the 
baselines  for  the  different  data  sets  was  as  follows:  GEM9 
( ±17cm) ,  GEM10B  (±llcm),  PGSL1  (±9cm.)  Rapp  (±9cm). 

The  solution  of  this  paper  performs  nearly  as  well  as  the 
PGSL1  solution.  This  is  due  to  the  inclusion  of  the  PGSL1  coef¬ 
ficient  set  into  our  solution  with  fairly  large  weights. 

The  next  test  was  carried  out  using  Starlette  Laser  data. 

Here  approximately  six  5  day  arcs  (with  2h  day  overlap)  were 
analyzed  using  the  Rapp  field  to  36,36  only.  The  following 
are  the  average  differences  using  the  specified  mode:  PGS1331 
( ±0 . 65m) ,  GEM10B  (±2.1m),  GEM9  (±2. 6m),  PGSL1  (±2. 5m),  and 
Rapp  (±4. 2m).  Clearly  the  solution  of  this  paper  doesnot  work 
as  well  as  the  other  solutions. 

Another  test  was  made  of  the  11  order  coefficients  using 
the  lumped  11  order  coefficients  given  by  Wagner  and  Lerch  1978). 
This  test  is  of  special  interest  as  no  specific  11  order  reson¬ 
ance  coefficients  were  incorporated  in  our  a  priori  potential 
coefficient  set.  The  root  mean  square  difference  was  computed 
for  14  separate  lumped  coefficients  as  given  by  Wagner  and  Lerch 
and  as  computed  from  several  coefficient  sets.  The  results 
are  given  in  Table  10. 


Table  10 


RMS  Lumped  Coefficient  Comparison 
(Wagner  and  Lerch  (1978))  minus  Computed  Value) 


Potential  Coefficient  Set 


GEM 9  (to  n=20 ) 
PGSL1  (to  n=20) 
PGS1331  (to  n=20 ' 
PGS1331  (to  n=36 ; 
Rapp  (to  n=20) 
Rapp  (to  n=36) 
GRIM3  (to  n=20 ) 
GRIM3  (to  n=36 ) 


RMS  Difference  (X109) 


±13.9 

±13.9 

±15.6 

±15.4 

±12.5 

±13.5 

±22.9 

±22.4 


And  our  final  comparison  was  with  a  set  of  24  hour  satellite 
accelerations  determined  by  Wagner  (1982,  private  communication). 
Using  8  different  observed  accelerations,  Wagner  tested  the 
following  coefficients:  (2,2),  (3,1),  (3,3),  (4,2),  (4,4), 

(5,1),  (5,3),  (5,5),  (6,2),  (6,4),  (6,6).  The  root  mean  square 
difference  between  the  observed  and  the  computed  acceleration 
(in  units  of  10~8rad/day  )  are:  GEM9  (c5),  GEM10B  (±3.6),  Rapp 
(±2.5) ,  and  PGSL1  (±1.0) . 

The  point  to  be  made  here  is  that  these  tests  indicate 
that  this  new  solution  for  some  cases  is  better  than  some  existing 
solutions.  Additional  testing  is  needed  to  obtain  a  more  complete 
picture  of  the  performance  of  this  new  coefficient  set  in  orbital 
work.  One  should  not  expect,  this  model  to  compete  with  models 
that  have  been  specifically  tailored  to  a  given  satellite. 


Summary  and  Conclusions 


We  have  generated  new  gravity  field  models  based  on  improved 
theory  and  improved  data.  We  have  used  more  current  satellite 
models,  terrestrial  data,  and  Seasat  altimeter  data.  The  data 
used  included  l°x  1°  data  so  that  the  sampling  error  could  be 
reduced.  We  have  incorporated  in  our  error  analyses  the  sampling 
error,  as  obtained  by  Colombo.  This  has  enabled  us  to  asses 
the  accuracies  of  the  various  components  of  the  model,  such 
as  undulations  and  anomalies. 

The  adjusted  anomalies  in  our  solutions  were  expanded  into 
spherical  harmonics  to  degree  180  ana  for  some  applications 
to  degree  300.  The  expected  error  in  the  coefficients  reaches 
almost  100%  at  the  higher  degrees  (120). 

The  adjusted  coefficients  were  used  to  compute  anomalies 
which  were  compared  to  the  input  l°x  1°  data  set  for  anomalies 
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whose  standard  deviation  was  ±7  mgals  or  smaller.  The  result 
for  the  adjusted  coefficients  and  several  other  coefficient 
sets  is  given  in  Table  11 . 


Table  11 

Mean  Square  Difference  Between 
Input  l°xlc>  Anomalies  and 
Anomalies  Computed  from  the 
Potential  Coefficient  Set 


Field 

NMAX  =  20 

NMAX  =  36 

GEM9 

210  mgal2 

•  — 

PSGL1 

210 

— 

PGS1331 

200 

183 

SET1 

202 

190 

Ad  j  . 

182  mgal 2 

180  mgal2 

We  see  a  slightly  better  agreement  for  the  adjusted  field, 
but  it  is  not  substantial.  The  reason  for  this  is  the  low 
relative  weight  assigned  to  the  l°x  1°  anomaly  data.  The 
adjusted  coefficients  to  degree  50  are  given  in  the  Appendix. 
The  complete  set  to  degree  180  (or  to  300)  is  available  on 
tape  as  are  the  adjusted  l°x  1°  anomalies. 

At  the  conclusion  of  this  study  it  is  clear  that 
several  things  could  have  been  done  differently.  Some  items: 

1.  Sea  surface  topography  corrections  should  be  made 
to  the  altimeter  data  to  reduce  long  wavelength  errors  in 
the  derived  anomaly  field. 

2.  Better  terrestrial  data  is  needed.  This  is  especia 
true  in  areas  where  geophysical  anomalies  now  exist. 

3.  The  optimization  procedures  for  the  combination 
of  these  data  types  should  be  implemented  to  assume  a  more 
rigorous  combination  solution. 
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Fully  Normalized  Potential  Coefficients 
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0.0029 

0.003d 

29 

lb 

O.OOCo  —0.0191 

0.002b 

0.0031 

29 

Id 

-0.0032  —0.0015 

0.0029 

0.0032 

29 

20 

— C  .0072  -U.0012 

0.0028 

0.0029 

29 

Zc 

C . 0129  0 • 0  Ol  3 

U. 003  2 

0.0021 

29 

29 

0.0010  —0.0  093 

0.0021 

0 .0026 

29 

2b 

0.0058  -0.0132 

0.0032 

0  .002  8 

29 

2d 

0 . outit  -0.0053 

0.0029 

0.0021 

30 

1 

0.0019  0.0090 

0.0021 

0  .001b 

30 

3 

-0.0009  — 0.0 l9fc 

0. 001 b 

0.0032 

30 

5 

0.0003  -0.0 Uq3 

0.0028 

0.0032 

30 

7 

0.0079  0.000c. 

0 .0030 

0  .0021 

30 

9 

-0.0073  -0.0039 

0.002b 

0.0019 

30 

11 

-0.00«>2  0.0  095 

U.G029 

0.0030 

30 

13 

0 .00d9  O.OOlo 

0.002b 

0.0021 

NOTE:  ALL  VALUES  TO 


CO 

19 

0.0021 

0.0065 

0.0032 

0.0031 

la 

lb 

0.0066 

-0.0110 

0.0039 

0.003o 

CO 

18 

-0.013b 

0.0078 

0.0037 

0.0034 

lb 

20 

0.0037 

-0.0127 

0.0031 

0.0035 

2b 

22 

0.0116 

O.OObb 

0.0035 

0.0031 

CO 

29 

0.0052 

0.0129 

0.003b 

0.0029 

CO 

2b 

-0.001U 

0.005b 

0.0023 

0.0030 

cl 

2 

0.0015 

0  .0002 

0.0019 

0.0034 

27 

9 

-0.0030 

0.0091 

0.0031 

0.0023 

2  7 

b 

0.0055 

0.0096 

0.0028 

0.0033 

27 

8 

0.0092 

-0.0196 

0.0029 

0.0032 

27 

10 

-0.0102 

0.0020 

0.0033 

0.0034 

cl 

12 

-0.0059 

-0.0062 

0.0030 

0.0029 

cl 

14 

0.0076 

0.0110 

0.0034 

0.0029 

2  1 

1  a 

0.0097 

-0.0038 

0.003b 

0.0021 

Cl 

18 

-0.0060 

0.0109 

0.0017 

0.002b 

27 

20 

-0.0003 

0  .0013 

0.0023 

0.0015 

27 

22 

-0.0099 

0.0017 

0.0022 

0.0032 

27 

2* 

-0.005o 

0.0007 

0.0025 

0.0015 

27 

2b 

-0.0091 

-0.00  59 

0.0029 

0.0030 

2d 

1 

-0.0016 

0.0052 

0.0012 

0.0026 

28 

3 

0.0002 

0.0091 

0.0009 

0.0031 

2b 

0 

0.0057 

-0.0110 

0.0022 

0.0034 

23 

7 

-0  .0007 

0.0029 

0.0033 

0.0016 

28 

9 

0.0082 

-0.0091 

0.0031 

0.0030 

*8 

11 

-0.0001 

-0.0017 

0.0015 

0.0034 

26 

13 

0.0050 

0.0096 

0.0034 

0.0025 

28 

15 

-0.0157 

-0.0051 

0.0032 

0.002* 

26 

17 

0.0111 

-0.002  7 

0.0034 

0.0031 

26 

19 

-0.0012 

0.0198 

0.0021 

0.0034 

26 

21 

0.0105 

0.0033 

0.0032 

0.0017 

t-6 

23 

0.003o 

-0.0025 

0.0029 

0.0033 

26 

25 

0.0010 

-O.Olbb 

0.0016 

0.0024 

c& 

27 

-0.0057 

0.00*5 

0.0032 

0.0022 

29 

1 

-0.0017 

-0.0031 

0.0031 

0.0029 

29 

d 

0.0012 

-0  .00  69 

0.0009 

0.0025 

c9 

5 

-0.0021 

C  .00 13 

0.0014 

0.0012 

29 

7 

-0.0016 

-0.00b3 

0.0025 

0.0032 

29 

9 

-0.0005 

0.003b 

0.0014 

0.0020 

29 

11 

— C .0037 

0.0057 

0.0025 

0.0032 

29 

13 

-0.0023 

-0.0083 

0.0030 

0.0012 

29 

15 

-0.0072 

-0.Q098 

0.001b 

0.0017 

29 

17 

0.0039 

-0.0O18 

0.0008 

0.002* 

29 

19 

-0 .OObb 

O.OObb 

0.0030 

0.0029 

29 

21 

-0.0039 

0.0002 

0.002b 

0.0016 

29 

23 

—0.007* 

0.0011 

0.0033 

0.0017 

29 

25 

0  .0092 

0.00b9 

0.0031 

0.0026 

29 

27 

-0.0067 

-0.0015 

0.0027 

0.0026 

29 

2? 

0.0073 

-0.0078 

0.0023 

0.002b 

30 

2 

-0.0131 

0.0007 

0.0032 

0.0014 

30 

9 

-0.0018 

0.0009 

0.0019 

0.0027 

30 

b 

0 . 003'. 

0.0061 

0.0027 

0.0025 

30 

6 

0.0000 

0  .0107 

0.001b 

0.0029 

dO 

10 

-0.0017 

-0.0096 

0.0021 

0.0019 

30 

12 

0.0150 

-O.OObO 

0.0032 

0.0025 

30 

19 

0  .0069 

0.00b9 

0.0030 

0.0020 

MULTIPLIED  3  Y  10**6 
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_ M_ 

CiL.M)  la  t  L  »M)  blOMA.  L 

S  I  bM  A  £ 

.  L 

_ iL 

_ L IL4M)  _  S (L*  M  )  SIGMA 

~J t' 

13* 

-0.U0.3  — vj.0006  0  .0014 

(J.003  1 

bO 

lo 

-3 toob 3 0.000/  oTouSv  DTooTS 

J.0 

17 

-C  . 003  0 

-0.0034 

0.GC2  1 

0  .002-* 

30 

io 

-0.0090 

-0.007? 

0.0030 

0.0024 

3u 

19 

— C  . GQ9o 

0.0013 

0 .002  7 

0  .0  0  i  6 

30 

40 

-0.0039 

0.0049 

0.0022 

0.0027 

3u 

21 

-0  .uU5d 

-0.0O30 

C .002  3 

0.0022 

60 

22 

-0.C013 

0.0005 

0.0013 

0.0020 

Ju 

23 

0.0027 

— b.OOob 

0.0024 

0  .0029 

JO 

24 

-o.uoio 

-0.003o 

0.0028 

0.0024 

30 

23 

0.0019 

-0 .0 1 3b 

0.C01 3 

0.0031 

jO 

2  o 

0.0012 

0.0106 

0.0032 

0.0030 

3o 

27 

-C.OOcl 

0.0  07  4 

0.0 02  7 

0.0014 

30 

2  S 

-0.0091 

-0.0135 

0.0016 

0.0017 

30 

29 

0.0030 

0.00**9 

G.0013 

0.0021 

30 

30 

-0.0033 

C  .007  5 

0.0010 

0.0010 

31 

1 

0. OGob 

—0 .0  Ob  3 

0.002b 

0.0043 

31 

2 

0.003? 

0.0029 

0.0024 

0.0015 

3i 

3 

*4J  .Cb3  1 

—0 . 0  Ob  o 

0.0017 

0.0027 

31 

4 

0.0066 

-0.0059 

0.0029 

0.0008 

31 

3 

-O.0013 

0 .0  022 

0.0024 

U  .0017 

31 

o 

-0.0016 

0.0025 

0.0016 

0.0010 

bl 

7 

0.0017 

-0.0014 

O.C024 

0.0010 

31 

b 

0.0011 

0.0002 

0.0003 

0.0016 

i  a 

9 

.bOl? 

O  .  0  04** 

0.0014 

0.0004 

3i 

10 

-0  .0029 

-C .00  37 

0.0008 

0 .0019 

31 

11 

-O .000 1 

0.0  lbV 

0.001  1 

0.003i 

3  1 

12 

0  .002b 

Q.00o5 

0.0023 

0.0029 

bl 

13 

0  .0 Gob 

0 .0037 

C . 0032 

0-002** 

31 

1** 

-0 .0050 

0.0035 

0.0028 

0 .0028 

31 

13 

0 .0063 

-o.uote 

0.0022 

0.0019 

31 

lo 

-0.0093 

0.005o 

0.0029 

0 .002  7 

31 

17 

-O  .  00  3b 

w  .00**3 

0.0020 

0  .0046 

3  i 

lb 

-0.0020 

-0.0011 

0.0026 

0.0020 

31 

19 

0.001b 

0  .o  ul  3 

C.C022 

0.0023 

31 

20 

-0.0009 

0.0031 

0.0008 

0.0030 

31 

21 

-O.0033 

C .  0  ii*t 

0.0024 

0.0014 

31 

22 

-0.0068 

-0.0041 

0.0018 

0.0017 

31 

23 

O.OOdl 

0.0073 

0.0023 

0.0024 

3  1 

24 

-0.0011 

-0  .OOOo 

0.0015 

O.OOlb 

31 

23 

-0  .0103 

-O.OulS 

0.002b 

0.0030 

31 

26 

-0.0093 

0.0006 

0.0023 

0.0028 

3i 

27 

O.002t 

0.0  13<* 

L .0030 

0.0032 

31 

23 

0.0079 

0.0022 

0.0029 

0.0029 

31 

29 

-0 .  oo3b 

-0.0043 

0.0013 

0  .0031 

31 

30 

-0.0016 

-0.003? 

0.0003 

0.0025 

31 

31 

-0.u0**d 

0.0 004 

0.0014 

0  .0002 

32 

1 

-0.0021 

0.0049 

0.002  6 

0.0019 

3^ 

2 

0  .0036 

-0.0033 

0.0023 

0.0023 

32 

3 

-0.0052 

0.0023 

0.0020 

0.0005 

32 

4 

0.0003 

-0.0091 

0.000b 

0  .0030 

32 

6 

0.0030 

0.0001 

0.0010 

0.000b 

32 

o 

-0 .0002 

-0 .0  033 

0.0004 

0  .0022 

32 

7 

-0.0023 

0.0003 

0.0023 

0.0004 

32 

a 

0.0093 

0.0072 

0.0030 

0  .0027 

32 

9 

0.0024 

-o.ooib 

0.0018 

0.0023 

32 

10 

0.OO14 

-0.00o2 

0.0023 

0  .0023 

64 

11 

-0.0025 

0.0037 

0.0012 

0.0021 

32 

12 

-0.0093 

0.0  l**fa 

0 .0029 

0  .0030 

32 

13 

0  .ooo> 

0.0075 

0.0012 

0.0026 

32 

l4* 

-0.0037 

0 . 0  C  12 

0.0027 

0.004** 

32 

16 

0.004b 

-0.0051 

0.002  3 

0.002b 

32 

lb 

0.O005 

-o».0oi2 

0.0024 

0  .0017 

32 

17 

-0.0053 

0.006b 

0.0G26 

0.0019 

32 

lb 

0  .o0b3 

0.0004 

0.0021 

0.0014 

32 

19 

-0.0005 

-0  .0006 

0.0006 

0.0024 

32 

20 

O.L0*»0 

0.0023 

0.0023 

0.0028 

32 

2  1 

-0.002? 

0.009b 

O.QClb 

0.0026 

32 

22 

—0 . bOd  3 

-u.0033 

0.0023 

0  .0029 

34 

23 

0.G04o 

-0.0018 

0.002b 

0.0016 

32 

2** 

—0 .0092 

0.0033 

0 .0026 

0.0021 

32 

23 

-0.0119 

-0.0003 

0.0026 

0.0026 

32 

2o 

0.0016 

—0 . 0  0 1 2 

O. 002  3 

0.0023 

22 

27 

-0.0065 

-0.0039 

O.OOlb 

0.0025 

32 

2d 

0.0032 

-0.0019 

0.0021 

0.0020 

32 

29 

-0.0010 

0.0041 

0.002  5 

0.0027 

32 

3L 

-o  .  006  7 

0.0 00 3 

o.002b 

0  .0017 

34 

31 

-0.0043 

0.0002 

0.0024 

0.0007 

32 

32 

O.CO**5 

0.C010 

0.001** 

0.0020 

33 

1 

-0.0029 

-0.0030 

0.002  3 

0.0026 

33 

2 

-0.0030 

0.0  Ol 3 

0.UO2O 

0  .001  i 

33 

3 

-0.0033 

0.0022 

0.0026 

0.0018 

33 

4 

—0 .0019 

0.0009 

0.0027 

0  .0016 

33 

3 

-0.0003 

0.0068 

0.0010 

0.0020 

33 

b 

-O.OOUO 

-0.0013 

C.COO** 

0.0004 

33 

7 

-0.0071 

-0.0002 

0.002  7 

0.0003 

33 

d 

0 .6012 

0  .Uu8*t 

U.UUG7 

0  .00  2o 

33 

9 

-0.001b 

0 .0020 

0.0019 

0.0014 

33 

10 

-0  •  0046 

-0.0033 

0.0027 

0  .0017 

33 

11 

-0 .0060 

-0.0062 

0.0018 

0.0027 

33 

12 

-c.uoio 

0.0064 

0.0010 

0.002b 

33 

13 

0.0043 

0.0034 

0.0027 

0.0012 

33 

14 

-t • ooO 1 

0.0030 

0.0012 

0  .0017 

33 

16 

-0.001b 

-0.0062 

0.0021 

0.0024 

33 

lb 

O.oOiO 

o.om 

0.0017 

0  .0008 

33 

17 

-0.0004 

0.0075 

0.0010 

0.0016 

33 

lo 

-U  .0043 

-0.0  03  1 

O.0C12 

0  .0022 

33 

19 

0.0066 

-0.0014 

0.0022 

0.0015 

33 

20 

0.0003 

—0 . 0  0  i  9 

O.OUlD 

0  -0013 

33 

21 

0.0011 

0.0056 

0.0007 

0 .002  7 

33 

22 

—0  •  CC**7 

-0.0 130 

0.0020 

0  .0028 

63 

23 

-0.0023 

-0.0110 

0.0015 

0.0028 

33 

24 

0.007o 

-0.0036 

0.002b 

0  .002  7 

33 

25 

0  .0034 

-0.0083 

0.0024 

0.0024 

33 

2b 

0.0077 

0 .0  Obb 

O.Ou27 

0  .0023 

33 

27 

-0.0031 

-0.0034 

0.0029 

0.0C24 

33 

2b 

-0  «OGO*t 

—0 .0  U  14 

0.0024 

0  .0023 

33 

29 

-C  .0173 

0.0065 

0.0019 

0.0031 

33 

30 

-0  .O01 6 

-0.0 149 

0.OC11 

0.0007 

33 

31 

0.0033 

-0.0003 

0.0015 

0.0002 

33 

32 

0 • 0G54 

-0.003b 

0.0028 

0  .0021 

33 

33 

0.0019 

0.0042 

0.0020 

0.0020 

NCjTc  :  ALL  VALUtS  TO  8t  MULTIPLIED  BY  10**6 
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_k. 

_ £3. 

--lik-iilL 

»mi  sigma  C  SIGMA  S 

_ L. 

_f± 

CtL.M) 

>(L.M) 

SIGMA  C  SIGMA  S 

34 

i 

— O.uOlb 

O.uuOb 

o.ooia 

0  .002b 

b4 

2 

C.LO04 

O.OOlb 

0.0G61  0 .061o 

34 

3 

0.0123 

O.LObb 

0.0026 

0  .0027 

b4 

4 

-0.002b 

-0.0002 

0.0015 

0.0005 

34 

b 

-0.U027 

0  .004.0 

0.0013 

0  .000^ 

34 

6 

0  .005© 

0  .0062 

0.0017 

0.0025 

34 

7 

-O  .UUOb 

—0  .001b 

O.OOli 

0.000b 

b  4 

6 

-0.011b 

0 .0022 

0.0029 

0.0011 

34 

9 

O.OOOh 

0.0008 

0.0003 

0  .0017 

b4 

10 

-0.006b 

-0.0012 

0.0026 

0.0020 

3* 

11 

-0  . O020 

-0.0  011 

C.COlo 

0  .001b 

b4 

12 

0.0101 

-0.0044 

0.0027 

0 .002© 

b4 

13 

-0  .0015 

C .0030 

0.002b 

C .0020 

34 

14 

-C.0031 

0.0096 

0.002b 

0.002© 

34 

lb 

C . 0002 

0. 0  Ot  2 

0  •  0  02  3 

0.0024 

b4 

1© 

-0.0024 

-O.OObb 

0.0009 

0.0018 

34 

17 

-0.0031 

U.OOOb 

0.0015 

0 .0007 

34 

16 

-0  .0097 

-0  .0040 

0.0023 

0.0019 

34 

19 

O.uOlb 

0 .0  029 

0.0023 

O.OOlo 

34 

20 

0  »C04b 

-0.0038 

0.0020 

0.0021 

3* 

21 

O.COGi 

-0.0052 

0.0012 

O.OOlb 

34 

22 

-0  .0018 

0.0052 

0.0021 

0.0019 

34 

23 

O.OuOb 

-0.0  006 

0. O004 

0  .002b 

34 

24 

0.007o 

0.0014 

0.0020 

0.0021 

3-4 

Cb 

O.005O 

-O.0077 

0.0017 

0  .0027 

34 

2b 

0.0023 

-C.0094 

0.0023 

0.0023 

34 

27 

0.0115 

-0  .0022 

0.002b 

0.0022 

34 

2  o 

-0.0001 

-0.0192 

0.0024 

0.0024 

34 

29 

0 • 00  bo 

— 0  .0ub3 

0 .002  o 

0.002b 

34 

30 

-0.0229 

-  0 .0006 

0.0028 

0.002b 

34 

31 

-0 . 00 1 1 

0  .  o  ©02 

0.0009 

0.0016 

34 

32 

0.0044 

-0.0002 

0.0024 

0.0002 

34 

33 

0.0120 

0.UU24 

0.0026 

0.0011 

34 

34 

-0 .0051 

-0.0004 

0.002  3 

0.000  8 

3b 

X 

-0.0032 

-0.Ui.lU 

0.0022 

0 ,0U2o 

bb 

2 

-0.0112 

0.0017 

0.0025 

0.0015 

3b 

0 .0060 

o.oooi 

U.0C2b 

0 .0014 

35 

4 

-0.001© 

0.0015 

0.0016 

0.0020 

3b 

b 

—0  *  0044 

-0.0  02  7 

0.0016 

0.0020 

3b 

O 

0.0032 

0  . 00©4 

0.0023 

0.002© 

3b 

7 

-0  .0022 

0.0022 

0.0022 

0 .0003 

b5 

a 

0  .0021 

0  .0026 

0.0006 

0.0010 

3b 

9 

-0.0004 

—0 .000  o 

o.OCll 

0 .00 16 

3b 

1C 

-0.0035 

0.0049 

0.0014 

0.0023 

3b 

11 

0.0010 

-0.0012 

0.0019 

0.0020 

35 

12 

0.0047 

-0.0041 

0.0022 

0.0026 

35 

13 

0.0020 

O.OUbo 

0.0018 

0.0021 

bb 

14 

-0.0054 

-0.0057 

0.0027 

0.0022 

35 

lb 

-0.O135 

O.OObb 

0.0024 

0.0019 

35 

1© 

-0.0044 

-0.0056 

0.0024 

0.0022 

3b 

17 

C.0022 

-0.0 068 

0.0011 

0 .0019 

3b 

16 

-0.0024 

-0.0058 

0.0006 

0.0025 

bb 

19 

-0.0004 

—0.0  02  7 

0.0018 

0 .0006 

b5 

20 

0  .0006 

0.00  16 

0.0014 

0.001© 

35 

21 

0. C085 

U.OOOb 

O.0027 

0.0012 

bb 

22 

0.0011 

0.0009 

0.0007 

0.0012 

3b 

23 

-0 .0043 

-O.0000 

0.0018 

0.0014 

b5 

24 

0  .0021 

0.0011 

0.001© 

0.0014 

3b 

25 

0.0062 

0.0Cb4 

0.0026 

0.0016 

35 

26 

-0.0041 

-0.0002 

0.0014 

0.0020 

3b 

27 

0 . oOoO 

—0.0 ibb 

0.0021 

0.0029 

35 

28 

0.0107 

-0.0160 

0.0028 

0.0029 

3b 

29 

0.0097 

0.002b 

0.0016 

0 .0023 

35 

30 

-0.0039 

0.00b9 

0.001© 

0.0021 

3b 

31 

0.0016 

0.0024 

0.0012 

0.0020 

bb 

32 

-0.0014 

-0.00©7 

O.OOlb 

0.0024 

35 

33 

0.0035 

— 0.0  ul9 

C.0010 

0.0012 

bb 

34 

-0.0018 

0.0018 

0.0010 

0.0001 

35 

3b 

—0  .00  b 6 

-0.0044 

u.0023 

0.0025 

bb 

1 

0.0024 

0.0033 

0.0024 

0.0023 

36 

2 

-O .00 (2 

-O  .OOue 

0.002  b 

0.0024 

36 

3 

0.0013 

-0.0094 

0.0003 

0.0024 

36 

4 

0.0017 

-0.0012 

0.0006 

0.0019 

36 

5 

-0.0045 

-0.0000 

0.0017 

0.0020 

io 

6 

0.0136 

-0.001b 

0.0026 

0.0016 

36 

7 

-0.0004 

0.0059 

0.0004 

0.0024 

36 

o 

0.0001 

—0.0  036 

0.0019 

0.0020 

bb 

9 

0.0010 

-0.0002 

0.0010 

0.0008 

36 

10 

-0.0002 

0.0 0b2 

U.OOOb 

0 .0024 

3o 

11 

0.0005 

0.0003 

0.0021 

0.0020 

36 

12 

-0.0003 

-0.0047 

0.0004 

0 .002 1 

3© 

13 

-0.0053 

0.0038 

0.0025 

0.001© 

36 

14 

-0.0059 

-0  .0  Oo2 

0.0020 

0.0017 

bb 

lb 

-0.000© 

0.0019 

O.OOlb 

0.0013 

36 

lb 

0.0006 

— 0.0018 

C . 0008 

0.0007 

bb 

17 

0 • 006 1 

-0.0076 

0.0019 

0.0024 

bo 

16 

0.0010 

0.0049 

0.0014 

0.0023 

3o 

19 

-0.0055 

-0.0013 

0.0023 

0.0010 

36 

20 

-0.0040 

0.0  020 

0.0021 

0.0006 

3© 

21 

0.003© 

-0.0030 

0.0020 

0.0014 

36 

22 

O.oOO? 

— 0.CC18 

0.0020 

0  .0017 

bb 

23 

0.000b 

-0.0009 

0.0007 

0.0003 

36 

24 

u.0014 

-U  .0007 

0 .000  b 

0  .0016 

3© 

Zb 

-0.0013 

0.012© 

O  .00 1 6 

0.002o 

36 

4.6 

0.0044 

o  ,Oo4b 

0.0026 

0 .0020 

3© 

27 

-0.0083 

0  .0086 

0.0025 

0.0025 

36 

2d 

0.0018 

—0 .0032 

0.0013 

0.0016 

b© 

29 

0.0044 

-0.0014 

0.0024 

0.0026 

3© 

30 

-0 .OObO 

0  .OOnb 

0.CC26 

0.0023 

bo 

31 

-0.0070 

-0.0050 

0.002© 

0.002© 

3b 

32 

0.0056 

0.0011 

O.OOlb 

0  .0011 

3© 

33 

0 .0029 

-0.0031 

0.0014 

0.0023 

36 

34 

— Q.OObO 

0.0026 

0.0024 

0 .00 18 

36 

3b 

0  .0007 

-0.0085 

0.0010 

0.0025 

3b 

3b 

0.0013 

-O.OOlo 

O.OOlo 

0 .0015 

37 

1 

0.0005 

-0.0020 

0.0026 

0.0028 

37 

2 

0.0011 

-U.0 il 3 

0.002b 

0  .0026 

37 

3 

0.0001 

0.0006 

0.0028 

0.0028 

37 

4 

0.0057 

—0 . 0  O 1 4 

0.0028 

0.0023 

b7 

b 

-0.0060 

0.0079 

0.0026 

0.0028 

37 

6 

-u.0025 

0.0  uo2 

0.0026 

0.0026 

37 

7 

0  .0064 

0  .0064 

0.0028 

0.0026 

NOTE:  ALL  VALUES  TO  be  MULTIPLIED  BY  10**6 
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SIGMA,  1  SIGMA  5 

j  7 

o 

-0 .0033 

— u.ooio  O.uOlb  0  .Ou2o 

3  l 

9 

0 . UUGb 

-o  .004?  o70o2b  o7oo2d 

37 

it 

- C  . GOGu 

0 . 0  C  j  c. 

0  .0028 

0  .002b 

J? 

1  1 

0.0031 

C. 002b 

0.0028 

0.0028 

37 

14 

0 . tO  1  c 

U.oOOl 

0 . 002b 

C .002o 

j7 

1  J 

-0 .0002 

-0.0097 

0.0028 

0.002b 

37 

14 

-0  .01-32 

-U.OO20 

0.002  b 

0 . 002s 

j  7 

1  3 

0 .00?  3 

-0.0007 

0.002b 

0.0028 

37 

ib 

0.0022 

0.0  1.36 

O.Go28 

0 .002b 

37 

17 

0.0052 

-C.0039 

0.0028 

0.0028 

37 

lo 

b  .  00  c.  9 

0 .u  012 

0  .  0  o2  6 

0  .0028 

37 

19 

-0.0072 

0.0029 

0.0028 

0.0028 

37 

30 

-U  .0069 

—0 .u 0 39 

0.0028 

0  .0028 

3? 

4  1 

0.0031 

-0.0009 

0.0028 

0.002a 

3? 

33 

0.CCO4 

0.0o29 

0.0028 

0 .0023 

37 

23 

0.0011 

0.0006 

0.0028 

0.0028 

37 

44 

-0.003c 

-0  .00  39 

0.002b 

0 .Ou2b 

37 

25 

0.0033 

-0.0019 

0.0026 

0.002b 

j7 

40 

-0  .  GuU  1 

0 .OObb 

0  .O02b 

0  .0028 

37 

47 

-0.002b 

0.0050 

0.0028 

0.0028 

j7 

38 

G.Gi33 

O.OO96 

0.0028 

0  .0026 

37 

29 

0  .0083 

G .0046 

0.0C28 

0.0028 

37 

3U 

-0 .00c>3 

0.0 190 

0  .0028 

0 .002b 

37 

3  1 

0.0045 

-0.0082 

0.0028 

0.0028 

37 

34 

-O.OUU1 

0.0019 

C  .002  8 

0.002b 

3? 

33 

-0.0002 

-C.02 10 

0.0C28 

0.0028 

j7 

34 

0 . UG5  0 

0.0031 

0 . 0  C2  b 

0 .0026 

37 

3b 

-0.0071 

-0.0077 

0.0026 

0.0028 

37 

JO 

-0.U031 

-0.0  o33 

0.002c 

0 .0026 

j7 

37 

0.00b5 

-0.0035 

0.0028 

0.0026 

jb 

i 

0  .  u0*»  j 

0.0023 

0.0027 

0  .002  7 

36 

0  •  CCbj 

0.0041 

0.002? 

0.0027 

3o 

3 

0 . 0093 

0.0 00 7 

0.CC27 

0  .0027 

3o 

9 

0  .002b 

-0  .U0U6 

0.002? 

0.0027 

3o 

3 

-0  .uGod 

0.0 b7  / 

0.0027 

0 .0027 

36 

0 

-0.0137 

C  .0037 

0.0027 

0.0027 

3a 

7 

—0  . tO 1 1 

-0.0  cod 

0.0027 

0  .0027 

3a 

6 

0.0038 

C.  00 19 

0.0C27 

0.0027 

3o 

9 

C .0070 

-O.OoaS 

0.0027 

0 .0027 

3o 

10 

-0.0U30 

-0.0046 

0.0027 

0.0027 

ja 

11 

0 .Gold 

U.UUbO 

U.0O27 

0 .0027 

3  6 

1 2 

0.0015 

-0.0010 

0.0027 

0.002? 

JO 

1J 

-O.OO09 

— 0 . 0  0  9  7 

0.002? 

0  .0027 

3b 

1  9 

-0  .OObb 

0.0041 

0 .002  t 

0.0027 

3o 

U 

0.000b 

-0.0 00b 

0.0027 

0  .0027 

30 

lb 

-0  .0076 

0.0106 

0.002? 

0.0027 

jb 

17 

U.O031 

0 .0  09  3 

0 .002  ? 

0  .0027 

36 

1  6 

0.0095 

-0.0033 

0.0027 

0.0027 

a3 

19 

0  .0008 

-0 .0009 

c  .002  7 

0  .0027 

30 

20 

0.0022 

-0.0026 

0.0027 

0.0027 

j6 

21 

-0.0001 

-O.OOij 

C.C027 

0.0027 

30 

22 

0  .0017 

0.0064 

0.0027 

0.0027 

5o 

43 

0 .00 i c 

0 .0  090 

0.0027 

0.0027 

38 

29 

-0.008b 

0.0022 

0.0027 

0.0027 

JO 

43 

0.0003 

0.0033 

0.0027 

0  .0027 

3e 

2b 

-0.0054 

0.0031 

0.0027 

0.0027 

3b 

27 

-0.0039 

0.0096 

0.0027 

0  .0027 

3c 

20 

—0  .OObb 

-0.0008 

0.002? 

0.0027 

JO 

29 

0  .OOol 

0.0037 

0.0027 

0  .0027 

36 

30 

0.0005 

0.0025 

0.0001 

0.0022 

j8 

31 

0.b023 

-0.0033 

0.002  7 

0  .002  7 

3o 

34 

0.0059 

0.001b 

0.0027 

0.002? 

3o 

3  j 

C  .O0s J 

O.G IjO 

0.GC2? 

0.0027 

30 

34 

-0.0C70 

0.0009 

0.0027 

0.0027 

36 

33 

0 . 00b  3 

0 .0036 

0.0027 

0  .0027 

33 

3b 

0.0016 

C  .0007 

0.0027 

0.0027 

3b 

37 

-0 .uO j9 

O.QOib 

0.0027 

0 .0027 

68 

38 

0.0050 

-0.0014 

0.0027 

0.0027 

39 

1 

0.0006 

U • 0  099 

0.0026 

0  .002b 

39 

2 

0.0035 

0.0077 

0.0026 

0.0026 

39 

3 

-0.0033 

0.0039 

0.002b 

0 .002o 

39 

4 

-0.0078 

-0.0087 

0.0026 

0.0026 

39 

3 

O.00b3 

0.0073 

0.002b 

0.0026 

39 

6 

-0.0048 

0.0009 

0.0026 

0.0026 

39 

7 

0  .0017 

-0 . 003c. 

U.0026 

0.002b 

39 

8 

0.0008 

0.0098 

0.0026 

0.0026 

39 

V 

0.0073 

0.0O72 

U.0U2b 

0 .002b 

39 

10 

-0.0027 

0.0050 

0.0026 

0.0026 

39 

11 

0  .Olio 

-0.0017 

0.0026 

0.002b 

39 

12 

-0.0030 

0.0129 

0.0026 

0.0026 

J9 

lJ 

-0 .00 j9 

-U.0O19 

0.002c 

0  .002o 

39 

14 

-0.0045 

0.001? 

0.0026 

0.0026 

39 

A3 

-0. 0090 

0.0007 

0.002b 

0.002b 

39 

lb 

0.0013 

-0.0019 

0.0026 

0.0026 

39 

17 

-0  .001b 

-O.0007 

O.OC26 

0.002b 

39 

16 

0.0032 

0.0004 

0.0026 

0.0026 

39 

*9 

0.0037 

0  .0099 

0.002b 

0 .002o 

39 

20 

-0.0006 

-0.0067 

0.0026 

0.002b 

39 

21 

-0.0033 

—0 .0  0 19 

0.0026 

0 .0026 

39 

22 

-0.0077 

-0.0016 

0.0026 

0.0026 

39 

33 

-V  «O09O 

0.0 020 

0.002b 

0.0026 

j9 

24 

-0  .0081 

0.0079 

0.0026 

0.0026 

39 

23 

-0 .0099 

—0  .0  02  2 

0.002b 

0  .002c 

39 

2b 

-C  .002  7 

0.0090 

0.0026 

0 .0026 

j9 

27 

-0 . ouao 

-0.0 039 

C.0C26 

0.002b 

39 

28 

-O.OCtO 

-0.0087 

0.0026 

0.002b 

39 

29 

-0 .0092 

— O.0u2b 

0.002  6 

0.G02o 

39 

20 

0.0036 

-0.0113 

0.0026 

0.0026 

39 

Jl 

0. 0099 

-0.0 117 

O.002b 

0  .002b 

39 

32 

0.0007 

0.0069 

0.0026 

0.002b 

39 

Ji 

—O'.  00b i 

O.0030 

0 . 002  b 

0.002b 

39 

3*t 

0.0027 

0.0034 

0.0026 

0.002b 

39 

j  3 

-0.O19J 

0.0099 

0 .0026 

0  .002b 

39 

3b 

0  .0024 

-0.0019 

0.0026 

0.0026 

3v 

j7 

o . oO 1 o 

— 0 .0  034 

0.002b 

•  J  .  0  0  2  b 

39 

30 

-0.0027 

0.0034 

0.0026 

0.0026 

J9 

39 

C.CGOO 

-0.0 OOC 

0 . 000b 

0.0003 

90 

1 

0.000b 

0.0C07 

0.0026 

0.0026 

40 

2 

— O .0093 

0  .  U  02  9 

0.0026 

0  .002b 

90 

3 

-0.0004 

0.0007 

0.0026 

0.0026 

-♦0 

9 

0.002C 

-0.0073 

0.002b 

0  .002b 

90 

3 

0.0135 

-0.0029 

0.0026 

0.0026 

NOTE : 

ALL  VALUES  TQ  6 E  MULTIPLIED  BY  10**6 
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o 

-6.U6X3 

OTuaJA  G'.UoIb 

07uu7& 

aO 

7 

-3TC$i6 

-rtot  13  0.0026  070076 

*»0 

6 

0 .003  A 

0 . 0  026 

0.0u26 

0  .0026 

aO 

9 

-0 .0002 

0.0013 

0.0026 

0.0026 

4U 

lu 

-0  . O0a6 

o.uobc 

0 . 0  0  2  6 

C  .0026 

AO 

1  1 

0  . uOOo 

-0.001A 

0.0026 

0.0026 

*»u 

lx 

U.L03V 

U.U  Ox  a 

0.0C26 

0 . 0  0  2  o 

aG 

13 

-0.0057 

-0.0021 

0.0026 

0.0026 

*♦0 

lA 

-u . GO  b2 

0  .00x3 

0.0026 

0  .0026 

aO 

13 

-0  .0076 

-0.0017 

0.0026 

0.0026 

aO 

1  O 

-0.0022 

— 0 . 0  Oa  p 

0.0026 

0.0026 

aO 

17 

0.0030 

-0.0027 

0.0026 

0.0026 

*»C 

x3 

u . ouo2 

0 .0  Ai  A 

0.00X6 

0.0026 

aO 

19 

-0.0018 

-0.0042 

C .0026 

0.0026 

40 

20 

-0  •  OU2a 

0.0066 

0.0026 

0  .0026 

aO 

21 

-0 .0046 

-0.0001 

0 .002 6 

0.0026 

hO 

22 

-0.01  bV 

—0 .0 lto 

0.0026 

0 . 002  e 

'  aO 

23 

-0  .0014 

-0.01a3 

0.0026 

Q.002o 

•*0 

2a 

0  .0023 

0.0  uau 

0.002c 

0  .  0  0  2  o 

aQ 

13 

0.0008 

-0.003A 

0.0C26 

0.0026 

*♦0 

2o 

0.003S 

-0.000  i 

0.0026 

0.002o 

aO 

27 

-0.001  A 

0.0002 

0. 002  o 

0.0026 

aO 

28 

0.U028 

0.0ua7 

0.0026 

0.0026 

aO 

19 

0.001<6 

-0.0016 

0.0026 

0.0026 

•»C 

30 

0.0012 

— 0.0  Olb 

0.0022 

0.0022 

aO 

31 

-0.0016 

-0.00 11 

0.0026 

0.002o 

aO 

32 

-O.OOU7 

-0.0029 

0.0026 

0  .0026 

AO 

b  3 

-0.003a 

-0.0031 

0.002o 

0.0026 

*.0 

b4 

0.0023 

0.0  028 

0 .002  o 

0.0026 

aO 

35 

0  .0117 

-0  .0053 

0.0026 

0.0026 

aO 

36 

0.0009 

0 • 0  03 i 

0.0026 

0  .0026 

aO 

37 

-0.0071 

0.0C12 

0.0026 

0.0026 

*♦0 

36 

-0  . 000 A 

0  .  0  0  7  a 

0.0026 

0.0026 

aO 

39 

0 .0044 

0.0021 

0.002a 

0 .0022 

AO 

AO 

0.UO01 

-0.0  002 

0.0007 

0.0009 

a  1 

1 

-0.0015 

-0.0068 

0.0025 

0.0025 

Al 

2 

u . 003a 

C.Cc^a 

0.0025 

0.0025 

Al 

3 

0.001a 

0.0044 

0.0025 

0 .002  5 

tl 

A 

-O.OOAl 

0.00X6 

0.0025 

0  .0023 

Ai 

5 

0.0055 

-0.0015 

0.0025 

0.0025 

A  1 

o 

O.O005 

o.OolO 

u .002  3 

0.0023 

Al 

7 

0  .0004 

0  .0017 

0.0025 

0.0025 

•♦l 

3 

-C.CC3A 

—0  .  o  oa  5 

0.002  5 

0.0025 

Ai 

9 

-0  .0066 

0 . Q0a9 

0.0025 

0.0025 

Al 

10 

0.003a 

-O.GOUX 

0.0023 

0  .0023 

Ai 

11 

0.0013 

-0.0060 

0.0025 

0.0025 

Ai 

12 

—0.000a 

o  .0013 

0.0023 

0.0025 

Al 

lb 

-0.003a 

0.0030 

0.0025 

0.0025 

A  1 

lA 

C  .  v.0a3 

0 . 0  00 1 

0.0023 

u  .0023 

A  i 

15 

-C .0007 

0.0008 

0.0025 

0.0025 

Al 

16 

-0.0009 

-0.0056 

0.0025 

0.0025 

Al 

17 

-0.002& 

0.0012 

0.0025 

0.0025 

Al 

18 

—0 . O000 

Q.0oa5 

0.0023 

0.0025 

Al 

19 

-0  .0046 

-0.001A 

0.0025 

0.0025 

Al 

20 

O.LOao 

-0.0  069 

0 .002  3 

0.0025 

Al 

21 

0.0000 

-0.0006 

0.0025 

0.0025 

Al 

22 

-0  .0093 

—0.0  Ox:  3 

0.0023 

0.0025 

Al 

23 

0.0008 

-0.0150 

0.0025 

0.0025 

Al 

24 

0 .0069 

o.ouio 

0.0025 

0 .0025 

Al 

25 

0  .0006 

0.0030 

0.0025 

0.0025 

Al 

it> 

0.0045 

-O.0o71 

0.0023 

0.0023 

Al 

27 

-0.0017 

0.002a 

0.0025 

0.0025 

Al 

28 

-0.0011 

-0.006a 

0.0025 

0.0025 

Ai 

29 

-0.0052 

0.0058 

0.0025 

0.0025 

Ax 

30 

0.UO2A 

—0  .0  oa5 

0.0023 

0.0025 

Al 

31 

0.0108 

0.0011 

0.002  5 

0.0025 

Al 

32 

-0.0023 

O.G070 

0.0025 

0.0025 

Al 

33 

-0.005b 

0.0098 

0.0025 

0.0025 

Al 

3A 

-0.0019 

0.0003 

0.0023 

0.0025 

Al 

b3 

-0.01AA 

0.0080 

0.0025 

0.0025 

Al 

36 

0.0007 

0  .0  C'09 

0.0025 

0.0025 

Al 

37 

0.0016 

-0.0118 

0.0025 

0.0025 

Al 

38 

-O.GibS 

u .0  0x3 

0.0025 

0.0023 

Ai 

39 

-0.0009 

0.0000 

0.0012 

0.0006 

Al 

aO 

0 .oOOO 

-0  1 0  00  i 

0.0002 

0.0005 

Ai 

41 

0 .0026 

0.0058 

0.0020 

0.0023 

a2 

1 

-0 . 00a3 

0.0015 

0.002a 

0.002A 

A2 

1 

-0.0038 

-0  .00  34 

0.0024 

0 .0024 

42 

3 

O.C01S 

0.008a 

0  .GC2  a 

0  .002a 

a2 

A 

0.0018 

0.0015 

0.002a 

0.002a 

42 

3 

-0 .00o9 

-0.0091 

0.002a 

0.002a 

a2 

6 

0.0027 

-0.0032 

0.0024 

0.002a 

a2 

7 

0 .002  7 

-0.0066 

0.002a 

0  .002 A 

a2 

8 

0.0029 

—0  .0024 

0.002a 

O.C02a 

a2 

V 

— O .001a 

0.0019 

0.002  a 

0.002a 

a2 

10 

O.OOaO 

0.0036 

0.002a 

0.002a 

a2 

11 

0  .OOOo 

-0.0013 

O.0u2a 

0  .002a 

a2 

12 

0.0036 

-0-0062 

0.002a 

0.002a 

a2 

13 

-0.0023 

0  .0  o<:2 

0.002a 

0.002a 

a2 

1A 

-0.0032 

0  .0056 

0.002a 

0 .0024 

A^ 

lb 

-0.0025 

0 . C  Oo9 

O.0C2a 

0  .002a 

A  2 

16 

0  .0068 

-0.0044 

0.002a 

0 .0024 

a2 

17 

-0  .0033 

—0 .0  OaI 

C.002a 

0 .00  2A 

A2 

la 

-0.0072 

0.0052 

0.002a 

0 .0024 

a2 

IS 

-0.0063 

-0.0012 

O.002a 

0 .002a 

Al 

20 

0  .  GO?  8 

-C .0003 

0.0024 

0.002a 

a2 

21 
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0.C021 

0.0021 

HO 

7 

-0.0025 

0.0001 

0.0021 

0.0021 

Hb 

8 

O.uUii 

O.O02 1 

0.0021 

0.0021 

Hd 

9 

-0.0039 

G.OOhO 

0.0021 

0.0021 

Hb 

10 

-0.000H 

0.0 037 

0.0021 

0.0021 

Hd 

11 

0  .0002 

0.0016 

0.0021 

0.0021 

HO 

12 

0.000b 

—0 .00x0 

0.0021 

0.0021 

Hb 

13 

0.002b 

0.0013 

0.0021 

0.0021 

**8 

1 H 

— U.OOIh 

O.OOlb 

0.0021 

0.0021 

hB 

15 

0.U019 

-0.0028 

0.0021 

0.0021 

ho 

lo 

O.OoOl 

0.0011 

C.0021 

0  .0021 

Hd 

17 

0.002b 

0.0019 

0.0021 

0.0021 

*♦6 

Id 

-U.0021 

0.002b 

0.CC21 

0.0021 

Hb 

19 

-0.0015 

0.0022 

0.0021 

0.0021 

h8 

20 

— o.uo22 

0.0  Oh 2 

0.0021 

0.0021 

h8 

21 

-0.0012 

-0.0027 

0.0021 

0.0021 

•*3 

42 

-U.0002 

0.0013 

0.0021 

0.0021 

Hd 

23 

-O.OOCb 

-0.0023 

0.0021 

0.0021 

Hb 

iH 

-O.uODl 

-0.0 OCb 

0.0021 

0  .0021 

Hb 

25 

-0.0020 

-0.0021 

0.0021 

0.0021 

He 

2b 

0.O00J 

-0.0  Obi 

0.0021 

0.002i 

nd 

27 

-0.0063 

0.0053 

0.0021 

0.0021 

•*a 

id 

0 .OOHH 

-0 . 0  Ob  3 

0.0021 

0  .002 1 

Hb 

29 

-O.OOOh 

-0.00H8 

0.0021 

0.0021 

Hd 

30 

-0.0  02  i 

-O.OOlH 

0.0021 

0.0021 

Hb 

31 

-0.0019 

-0.0021 

0.0021 

0.0021 

HO 

3c 

0 .uOJo 

-0  .0023 

0.0O21 

0.0021 

HO 

33 

0.001S 

-0.0021 

0.0021 

0.0021 

Hd 

3h 

O.uOoh 

0  .OOH9 

0.0021 

0.0021 

Hd 

3  5 

-0.003b 

-0.0020 

0.0021 

0.0021 

Ho 

3b 

— U.OOlH 

0.0029 

0. 002 1 

0.0021 

Hb 

37 

-0.0031 

-0.0025 

0.0021 

0.0021 

H8 

3b 

-0.0101 

-0.0022 

0.0021 

0.0021 

Hd 

39 

0.0043 

-0.0090 

0.0021 

0.0021 

h3 

HO 

-C  .0001 

0.0 C3H 

0.0021 

0.0021 

Hd 

HI 

-0.0009 

-0.0095 

0.0021 

0.0021 

Hd 

h2 

0.0023 

0 . 0  03  3 

0.0021 

0.0021 

Hb 

h3 

0.0001 

0.0C58 

0.0007 

0.0020 

HO 

HH 

0.0023 

-O.OOOi 

0.0O21 

0.0021 

H3 

45 

0.0033 

0.00b7 

0.0021 

0.0021 

Hb 

Hb 

-0.0015 

o.OOdb 

0.0021 

0.0021 

H8 

H7 

0  .OObO 

0.0063 

0.0021 

0.0021 

Ho 

HO 

0.0031 

-0.0 OOd 

0.0021 

0  .0021 

h9 

1 

0.0033 

0.00C1 

0.0021 

0.0021 

H9 

2 

0.0U2H 

0*0032 

0.0021 

0.0021 

H  9 

3 

-0.0027 

O.OOHO 

0.0021 

0.0021 

H9 

H 

-U.0012 

0.0092 

0.0021 

0.0021 

H9 

3 

0.0010 

0.0004 

0.0021 

0.0021 

49 

b 

0 .0003 

0 . O  O2O 

0. 0021 

0  .0021 

h9 

7 

-0.0003 

0.0002 

0.0021 

0.0021 

49 

0 

-O.OOOi 

0.0023 

0.0021 

0.0021 

h9 

9 

-0.0002 

0  .00  70 

0.0021 

0.0021 

49 

10 

—0 • OObH 

C.OOOb 

0.0021 

0.0021 

H9 

11 

0.00H0 

0.0007 

0.0021 

0.0021 

tY 

12 

—0 • cOho 

—0  .0  0i3 

L.0021 

0.0021 

h9 

13 

0.0059 

0  .  OOhh 

0.0021 

0.0021 

49 

lH 

-o  .uooi 

-0.0  001 

O.OO2I 

0.0021 

H9 

13 

O.OOlb 

0.0008 

0.0021 

0.0021 

49 

lb 

-C .0012 

— 0.0  07 C 

0.0021 

0.0021 

H9 

17 

-O.OOlH 

-0.0018 

0.0021 

0.0021 

*»y 

io 

-O.OOlb 

-0.  OOOb 

0. 002 X 

0.0021 

h9 

19 

-0.0011 

-0.0021 

0.0021 

0.0021 

h9 

20 

0.0052 

0 .0  OUb 

O.0021 

0  .0021 

H9 

21 

-0.U012 

-0.0058 

0.0021 

0.0021 

h9 

22 

-O.001H 

O.0C2H 

0.0021 

0.0021 

H9 

23 

0.0050 

-0  .0007 

0.0021 

0.0021 

h9 

2h 

0 .00**3 

—0.0  oC  x 

0.0021 

0.0021 

H9 

25 

-0.0022 

0.002H 

0.0021 

0.0021 

H9 

2b 

-0 .0073 

-0.000  7 

C.C021 

0.0021 

h9 

27 

-0  .0027 

0.0022 

0.0021 

0.0021 

HH 

2o 

O.COld 

-u.0107 

C.CC21 

0.0021 

H9 

29 

-0.0020 

0.0033 

0.0021 

0.0021 

h9 

30 

0 • oOt . 

0.  0  01 1 

L.0021 

0.0021 

h9 

31 

U.0000 

-O.OGh? 

0.0021 

0.0021 

49 

32 

-0.0020 

-O.C031 

0.0021 

0.0021 

h9 

33 

0.0029 

-0.0035 

0.0021 

0.0021 

49 

34 

0.G0H2 

O.OOlH 

C.0021 

0.0021 

H9 

33 

0.0051 

0.0035 

0.0021 

0.0021 

hV 

3b 

-0 .  00h3 

0. 0  039 

0.0021 

0.0021 

h9 

37 

-0.0020 

0.0021 

0.0021 

0.0021 

49 

3d 

0 . 0039 

-O.OOlb 

0.0021 

0  .0021 

H9 

39 

0 .0040 

0.0023 

0.0021 

0.0021 

49 

hO 

—0 .002  0 

0.0  uGb 

0.C021 

0.0021 

h9 

HI 

-0.0017 

-0.0019 

0.0021 

0.0021 

h9 

h2 

-O. GQ3h 

b.OObS 

0.0021 

0.0021 

h9 

H3 

0.0050 

-C  .0097 

0.0021 

0.0021 

h9 

HH 

0 .00b 3 

0.007 1 

0.0021 

0.0021 

HV 

H3 

0  .0020 

-0.0002 

0.0021 

0.0021 

h9 

Hb 

0.0019 

O.UOHl 

0.0021 

0.0021 

h9 

H7 

-O.OOOh 

O.OOOI 

0.0021 

0.0021 

h9 

Hb 

-C.0OO3 

0.0017 

0 . C02 1 

0.0021 

H9 

H9 

0.0009 

O.OOhO 

0.0021 

0.0021 

3u 

1 

0.0013 

-O.002b 

0.0020 

0.0020 

30 

2 

-0.00b2 

-0 .00  bH 

0.0020 

0.0020 

30 

3 

-0.0001 

-0.0001 

C .0020 

0  .0020 

50 

H 

-0.0119 

0.002b 

0.0020 

0.0020 

50 

5 

-0.0017 

G.COCh 

0.0020 

0.0020 

30 

6 

O.OOOh 

0.0011 

0.0020 

0.0020 

50 

7 

0.C031 

O.OOH3 

0.0020 

0.0020 

30 

d 

-0.00h9 

-0.0026 

0.0020 

0.0020 

NOTE:  ALL  values  to  BE  MULTIPLIED  BY  1 0**6 


n8“=i&6l] 


5L 

11 

-0.0C3U 

C  .0  037 

C.0C20 

0 .0020 

50 

13 

G.uull 

0.0019 

0.0020 

0.0020 

50 

13 

-U.Q007 

-O.0035 

0.0020 

0.0020 

50 

17 

0  «u029 

-0.0  05  0 

C.CC2C 

0.0020 

50 

19 

u.0013 

O.OOifc 

O.0C20 

0 .0020 

30 

21 

-O.OOOo 

0.0  003 

0.0020 

0.0020 

50 

23 

-C.0Q16 

-0.0 Co9 

0.0020 

U.0020 

50 

25 

0.0063 

0  .0039 

0.0020 

0  .0020 

50 

27 

0  .0006 

-0.001*: 

0.0020 

0.0020 

50 

29 

0.0035 

O.Ottl 

L .0020 

0 .002  u 

30 

31 

0.0010 

0.0039 

0.0020 

0.0020 

50 

33 

— 0 .002** 

-0.003b 

0.0020 

0.0020 

5  C 

33 

0.t*02l 

0.0  0^3 

C.0020 

0.0020 

50 

37 

-L.O055 

0  .0003 

0.0020 

0 .0020 

50 

39 

—0.0035 

0.007V 

C.0020 

0.0020 

50 

tl 

—0 .00*8 

—0 .0  C**3 

0.002C 

0.0020 

50 

*3 

-0.0019 

-O.0046 

0.0020 

0.0020 

50 

45 

0.0000 

0.0032 

0.0020 

0  .0020 

50 

*»7 

-0.0U53 

—0.0  096 

0.0020 

0  .0020 

30 

49 

0.0031 

-0.0039 

0.0020 

0.0020 

NOTE:  ALL  V ALL) LS  TQ 


50 

12 

-0.0036 

0.00*0 

0.0020 

0.0020 

30 

14 

-0.0027 

0.0027 

0.0020 

0.0020 

30 

le> 

0  .0005 

-C  .0067 

0.0020 

0.0020 

50 

lb 

0.0031 

-0.0051 

0.0020 

0.0020 

30 

20 

0.0036 

-0.0004 

0.0020 

0.0020 

50 

22 

0.0006 

-0.0012 

0.0020 

0.0020 

30 

24 

0.0101 

-0.0000 

0.0020 

0.0020 

30 

2o 

-0.005b 

0.0022 

0.0020 

0.0020 

50 

28 

-0.0012 

0.0060 

0.0020 

0.0020 

30 

30 

0.0038 

0.0046 

0.0020 

0.0020 

30 

32 

-O.OOle 

0.0014 

0.0020 

0.0020 

30 

34 

0.0010 

-0.0009 

0.0020 

0.0020 

30 

36 

-0.0005 

0.0012 

0.0020 

0.0020 

30 

3b 

-0.0023 

-0.0096 

0.0020 

0.0020 

30 

40 

0.0043 

0.0055 

0.0020 

0.0020 

30 

42 

0.0032 

-0.0016 

0.0020 

0.0020 

30 

44 

-O.OOOo 

-0.0049 

0.0020 

0.0020 

30 

46 

-0.0039 

0.0061 

0.0020 

0.0020 

50 

40 

-0.0022 

0.0016 

0.0020 

0.0020 

50 

50 

0.0029 

0.0032 

0.0020 

0.0020 

&E  MULTIPLIED  BY  10**6 


